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ABSTRACT 


This  thesis  investigates  the  use  of  porous  materials  in  a  multi-layered 
armor  concept.  The  prototype  layered  structure  consists  of  an  initial  high- 
strength  material  to  slow  down  the  projectile  and  cause  significant  plastic 
deformation,  followed  by  an  orthotropic  wave-spreading  layer  to  spread  shock 
waves  laterally  away  from  the  axis  of  penetration  and  subsequently  attenuate  the 
shock  waves  by  using  a  porous  material  to  convert  kinetic  energy  into  internal 
energy.  Based  on  the  above  armor  concept,  composite  plates  consisted  of  an 
alumina  (A^Osj-based  ceramic,  Dyneema®  and  porous  foams  were  constructed 
and  tested  against  conventional  armor  steel  of  equivalent  areal  density.  This 
study  used  two  commercially  available  porous  materials,  one  based  on  aluminum 
metal  and  one  a  rigid  polyurethane  foam.  This  study  also  investigated  the  effect 
of  porous  initial  density  of  the  polymeric  foam  on  ballistic.  The  author  developed 
a  P-a  compaction  model  for  the  chosen  porous  materials  for  use  in  AUTODYN® 
simulations  to  describe  their  dynamic  compaction  during  an  impact  event.  The 
author  also  conducted  a  ballistic  trial  to  validate  the  performance  of  the  armor 
laminate  against  numerical  simulations.  Based  on  the  results  of  this  study,  the 
porous  layer  has  proven  to  be  a  good  shock  attenuator.  The  porous  material 
efficiently  delays  the  shock  wave  propagation  and  attenuates  the  amplitude  by 
absorbing  the  kinetic  energy  through  compaction  of  the  material.  The  current 
research  has  also  proven  that  the  material  layering  sequence  is  fundamentally 
correct  and  has  its  merits. 
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I.  INTRODUCTION 


A.  MOTIVATION 

With  ever-increasing  ballistic  threats  and  survivability  requirements, 
coupled  with  the  need  for  high  tactical  mobility  of  future  combat  systems, 
research  on  lightweight  armor  systems  continues  to  be  of  interest  to  many 
military  forces  around  the  world.  In  order  to  optimize  the  weight/performance 
ratio  of  armor  systems  for  both  personnel  and  vehicles,  it  is  important  to 
understand  at  a  very  fundamental  level  how  they  respond  to  ballistic  impacts. 
Achieving  a  reduction  in  armor  weight  coupled  with  improved  ballistic 
performance  will  require  the  use  of  new  materials  or  combinations  of  materials,  a 
better  fundamental  understanding  of  stress  wave  propagation  across  dissimilar 
materials,  the  optimization  of  competing  performance  requirements,  and  an 
improved  understanding  of  damage  mechanisms  of  new  ballistic  threats.  Taken 
together,  this  means  that  there  is  still  a  need  for  innovative  armor  concepts. 

There  has  been  a  considerable  evolution  in  armor  technology  for 
personnel  and  vehicles  since  World  War  II.  However,  most  passive  armor 
systems  have  focused  on  the  use  of  very  strong  materials  to  disrupt  the  incoming 
threat  by  breaking  up  or  deforming  the  projectile.  Hardened  steel  or  rolled 
homogeneous  armor  (RHA)  continue  to  be  popular  choices,  especially  for 
armored  vehicles  because  of  overall  good  ballistic  performance  and  its  relatively 
low  cost.  The  use  of  RHA  as  a  benchmark  for  new  armor  systems  still  occurs 
today.  However,  at  high  impact  velocities,  due  to  the  short  timeframe  and  highly 
localized  nature  of  the  impact,  the  material  strength  of  such  strong  materials  is 
vey  small  when  compared  to  the  compressive  and  tensile  stresses  generated 
upon  impact.  Hence,  there  can  be  considerable  damage  inflicted  on  the  target. 
To  effectively  defeat  a  projectile  without  experiencing  high  material  stress,  armor 
materials  designed  to  attenuate  impact  generated  shock  waves,  and  absorbs  the 
kinetic  energy  of  the  projectile  and  converts  it  to  another  form  of  energy,  i.e., 
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heat,  can  effectively  mitigate  the  penetration  process.  Porous  materials  possess 
good  shock  absorption  properties,  which  make  them  light  energy  absorbers, 
potentially  suitable  for  lightweight  armor  applications. 


B.  RESEARCH  GOAL 

The  goal  of  this  research  is  to  investigate  the  use  of  porous  materials  in  a 
layered  structure  for  advanced  armor  protection.  The  current  work  is  a 
continuation  of  thesis  work  completed  by  Poh  [1].  As  there  is  currently  no  single 
material  that  can  be  optimized  to  both  disrupt  and  absorb  the  kinetic  energy  of  an 
incoming  threat,  Poh  [1]  had  introduced  the  concept  of  an  integrated  multi¬ 
layered  structure  to  enable  sharing  of  functions  among  the  various  material 
constituents  to  optimize  ballistic  performance  for  improved  survivability.  The 
layered  structure  (see  Figure  1)  consists  of  an  initial  high-strength  material  to 
slow  down  the  projectile  and  cause  significant  plastic  deformation,  followed  by  an 
orthotropic  wave-spreading  layer  to  spread  the  shock  waves  laterally  away  from 
the  axis  of  penetration  and  subsequently  attenuate  them  by  absorbing  the  kinetic 
energy  through  compaction  of  the  porous  layer.  Poh  [1]  had  clearly  shown 
through  numerical  simulation  using  theoretically  desirable  material  properties  that 
the  layered  structure  concept  has  its  merits. 


Pene 


n  High  Strength  Layer 


Wave-spreading  Layer 


Porous  Layer 


Figure  1 .  Multi-layered  Structure  Armor  Concept 
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The  goal  of  this  research  is  to  better  understand  the  basic  shock  physics 
response  of  a  similar  concept  and  develop  a  hydrocode  model,  which  can  predict 
the  ballistic  performance  to  identify  material  properties  in  order  to  optimize  the 
armor  design  for  performance/weight  ratio.  In  particular,  this  study  can  address 
the  effect  of  pore  size  and  void  volume  of  the  porous  material  on  overall  ballistic 
performance.  In  future  research,  with  better  understanding  of  shock  wave 
propagation  across  dissimilar  materials,  there  could  be  the  development  of  a 
single  material  with  multi-functional  roles  to  integrate  the  roles  currently  fulfilled 
by  each  individual  material  layer. 

C.  RESEARCH  OBJECTIVES 

Since  the  work  done  by  Poh  [1]  has  shown  through  numerical  modeling 
that  the  multi-layered  structure  armor  concept  has  its  merits,  the  current  research 
objective  is  to  refine  the  design  and  prove  the  concept  experimentally  through 
ballistic  trials  using  commercially  available  materials  to  form  the  layered 
structure.  The  focus  of  this  thesis  is  on  the  porous  layer.  Concurrent  thesis  work 
by  Ong  [2]  focuses  on  the  wave-spreading  layer. 

The  author  selected  two  commercially  available  porous  materials,  namely 
aluminum  metal  foam  and  rigid  polyurethane  foam,  and  performed  material 
characterization.  He  also  investigated  Polyurethane  foam  of  two  different 
densities.  The  author  performed  numerical  modeling  using  the  AUTODYN® 
hydrocode  in  conjunction  with  ballistic  trials  to  better  understand  the  two- 
dimensional  response  and  shock  physics  of  the  porous  material  in  the  laminate 
configuration  described  above.  He  then  compared  and  validated  the  simulation 
predictions  against  the  field  trial  results  in  terms  of  deformation,  damage,  and 
ballistic  performance.  As  there  is  currently  a  wide  range  of  ballistic  threats,  the 
author  recognizes  the  limits  to  the  effectiveness  of  his  armor  design  against  each 
individual  threat.  Hence,  the  current  study  focuses  on  small  arms  threats  for 
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personnel  and  lightweight  vehicle  armor.  Then,  there  can  be  a  suitable  scaling  of 
the  proven  concept  to  meet  the  ever-changing  battlefield  survivability 
requirements. 


4 


II.  BACKGROUND 


A.  INTRODUCTION 

The  use  of  porous  materials  in  shock  absorption  and  isolation  applications 
is  extensive  because  of  their  effectiveness  in  attenuating  shock  and  mitigating 
impact  pressure.  A  porous  medium  or  material  is  a  solid  permeated  by  an 
interconnected  network  of  pores  (voids)  filled  with  a  fluid  (liquid  or  gas)  [3].  The 
pores  distend  the  solid  material,  causing  it  to  have  an  initial  density  lower  than 
the  crystalline  density.  Due  to  their  effectiveness  in  energy  absorption  for 
protection  applications,  such  materials  continue  to  be  of  interest.  However,  most 
current  applications  have  been  restricted  to  relatively  low  velocity  impacts. 

Recent  experimental  studies  [4]  have  demonstrated  that  a  porous  material 
such  as  polyurethane  foam  can  cause  peak  overpressure  attenuation  of  a  blast 
wave  by  10-50%  as  compared  to  air  as  a  medium.  Current  potential  military 
applications  include  using  rigid  polyurethane  foam  (RPF)  for  mine  protection. 
The  results  of  explosive  experiments  show  that  RPF  is  effective  in  absorbing 
direct  blast  waves  from  explosives.  Beyond  that  application,  the  foam  is  also 
effective  in  reducing  the  signature  of  vehicles  passing  over  anti-tank  (AT)  mines 
to  prevent  the  mine  from  firing.  Further  investigations  for  field  applications  are 
currently  on-going  [5].  Recent  experimental  work  on  low-density  aluminum  foam, 
as  a  backing  plate  for  ceramic  in  composite  armor,  has  revealed  improved 
ballistic  performance  with  the  foam  exhibiting  significant  non-linear  deformation 
and  stress  wave  attenuation.  The  superior  performance  and  unique  attributes  of 
porous  materials  may  allow  a  significant  advancement  towards  lighter  and  more 
damage-tolerant  composite  integral  armor  for  the  next  generation  of  lightweight 
armored  vehicles  [6]. 
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B.  DYNAMIC  COMPACTION  OF  POROUS  MATERIALS 


In  order  to  exploit  the  unique  properties  of  porous  materials,  it  is  important 
to  understand  their  behavior  when  subjected  to  dynamic  compaction.  Porous 
materials  respond  quite  differently  from  their  solid  parent  material  when 
subjected  to  a  shock  or  dynamic  compression.  For  solid  materials  (see  Figure  2) 
subjected  to  small  strains  and  stresses,  the  shock  Hugoniot  and  the  ambient 
isentrope  are  similar.  Hence,  for  qualitative  purposes,  one  can  assume  that  the 
Hugoniot  is  the  release  path  for  a  material  subjected  to  shock  compression  with 
only  a  small  error.  When  a  solid  material  is  subjected  to  a  shock  from  its  initial 
rest  condition  Po,Vo  during  the  shock  compression  process  from  initial  state  to  the 
final  state  P,V,  the  material  shock  response  follows  the  Rayleigh  line.  However, 
when  the  pressure  is  released,  the  material  follows  the  shock  Hugoniot  curve, 
which  is  an  approximation  following  the  released  isentrope,  down  in  stress.  The 
resulting  observation  is  that  the  shock  process  is  irreversible  as  the  internal 
energy  increase  due  to  the  shock  process,  which  is  proportional  to  the  area 
under  the  Rayleigh  line,  is  not  completely  recovered  upon  release  and  hence,  the 
material  is  left  with  some  residual  energy,  usually  in  the  form  of  heat  [7]. 


Figure  2. 


Shock  path 
Release  path 


Shock  Compression  Behavior  of  Solid  Materials 
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For  porous  materials  (see  Figure  3),  since  the  material  is  distended,  its 
initial  volume  is  larger  than  the  full  density  initial  volume.  The  inclusion  of 
porosity  in  the  matrix  of  the  solid  material  allows  for  the  attainment  of  larger  initial 
volumes.  When  subjected  to  shock  compression,  the  material  (assumed  here  to 
have  no  strength)  will  follow  the  Rayleigh  line  from  the  distended  initial  volume  to 
the  final  state  (diagonal  path),  and  if  the  shock  stress  is  high  enough,  approach 
the  full  density  solid  Flugoniot.  During  the  compression  from  1  to  2,  the  pores  in 
the  distended  material  are  completely  crushed  out,  and  this  path  is  known  as  the 
crush  curve.  The  pores  will  remain  crushed  once  crushed.  Flence,  when  the 
pressure  is  released,  the  material  will  approximately  follow  the  Hugoniot  down  to 
final  state  4,  which  has  a  smaller  ambient  volume  than  the  initial  distended 
volume  at  1  [7].  The  shock  compression  process  therefore  crushes  out  the  pores 
and  leads  to  permanent  material  densification.  The  release  path  from  state  2  is 
really  an  isentrope  centered  on  the  entropy  increase  that  happens  in  the  shock 
compression  process.  This  isentrope  will  be  close  to  the  Hugoniot  only  when  a 
small  entropy  increase  occurs  in  the  shock  compression  process.  The  author 
makes  this  assumption  here  to  be  able  to  qualitatively  understand  the 
shock/release  process. 
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Due  to  their  distended  volume,  as  compared  to  fully  dense  solid  materials, 
there  is  a  much  greater  internal  energy  increase  behind  a  shock  in  a  porous 
material  compared  to  a  solid  material  (see  Figure  4).  Porous  materials  can 
absorb  kinetic  energy  when  the  pores  collapse  so  their  internal  energy,  density, 
and  temperature  all  increase  through  p-v  work.  In  the  process  of  pore  collapse, 
the  material  can  strongly  attenuate  shock  waves,  and  the  kinetic  energy  is 
absorbed  and  turned  into  waste  heat  [7],  As  mentioned  earlier,  the  shock 
process  is  irreversible  and  some  fraction  of  the  internal  energy  increase  realized 
in  the  compression  is  not  recovered  upon  release  of  pressure.  Some  of  this 
energy  is  left  behind  as  waste  heat,  causing  the  final  temperature  to  be  greater 
than  the  initial  temperature.  It  is  a  well-known  fact  that  porous  materials  increase 
in  temperature  very  quickly  when  shocked,  much  more  quickly  than  solid  matrix 
material  alone.  This  can  cause  melting  and  even  vaporization  at  much  lower 
stress  than  for  solid  materials.  Hence,  it  is  clear  that  there  is  much  more  energy 
absorbed  or  residual  energy  left  behind  for  the  porous  material  compared  to  the 
full  density  solid  material.  The  presence  of  pores  makes  porous  materials  good 
mechanical  energy  absorbers. 


8 


A  stress-strain  diagram  illustrates  the  energy  absorption  characteristics  of 
porous  materials  [3].  The  stress-strain  plot  can  be  divided  into  three  regions  (see 
Figure  5): 


Eng.  Stress 


Region  1 

Region  1  is  known  as  the  linear  elastic  region.  Under  low  pressure  (or 
stress)  conditions,  only  elastic  deformation  occurs.  The  temporary  shape  change 
is  recoverable  after  the  load  is  removed  and  the  material  returns  to  its  original 
shape.  This  type  of  deformation  involves  only  stretching  or  elastic  compression 
of  the  atom-atom  bonds.  The  slope  of  the  plot  at  this  region  is  defined  as  the 
Young’s  modulus,  which  is  a  measure  of  the  stiffness  of  an  isotropic  elastic 
material. 

Region  2 

When  there  is  sufficient  stress  to  permanently  deform  the  material,  plastic 
deformation  occurs.  This  is  when  the  first  cell  wall  (or  pore)  collapse  occurs. 
Hence,  this  is  known  as  the  collapse  region.  When  the  material  is  compressed 
beyond  this  yield  point,  upon  release  of  the  applied  load  the  material  returns  to 
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zero  stress  (relaxed  state)  along  an  elastic  path,  resulting  in  a  residual  strain  or 
permanent  shape  change.  This  means  that  there  is  hysteresis  in  the  stress 
strain  response. 

Region  3 

As  the  porous  material  progressively  collapses,  densification  occurs,  i.e., 
volume  is  reduced.  This  is  known  as  the  densification  region.  Densification 
occurs  up  to  the  densification  strain  cd.  At  strain  values  beyond  eo,  the  structure 
is  completely  compacted,  the  stress  rises  steeply,  and  the  material  behaves  like 
a  fully  dense  material.  The  area  under  the  stress-strain  curve  (up  to  the 
densification  strain)  is  proportional  to  the  energy  the  porous  material  can  absorb 
and  is  the  energy-absorbing  performance  indicator  of  the  porous  material. 
Hence,  it  is  desirable  for  the  material  to  posses  a  large  eo  for  energy  absorption 
applications. 

Some  behavioral  properties  of  porous  materials  during  dynamic 
compaction  can  also  be  determined  analytically  using  conservation  equations. 
The  shock  and  particle  velocity  for  a  porous  material  can  be  determined  by  using 
the  jump  conditions  to  derive  a  relationship  between  pressure  and  run  distance 
or  time.  The  derived  relationship  reveals  that  the  pressure  drops  with  distance 
as  1/x^,  and  with  time  as  1/t.  This  shows  the  dispersive  nature  of  porous 
materials,  which  makes  them  good  shock  wave  attenuators  [7]. 

C.  COMPACTION  MODELS 

To  better  understand  the  behavior  of  a  shock-compressed  porous  material 
through  numerical  modeling,  constitutive  models  have  been  developed  to  model 
the  compression  of  porous  materials.  The  earliest  and  simplest  model  is  the 
snowplow  model  (relating  porous  material  compaction  to  snow  compaction).  In 
this  model,  the  assumption  is  that  there  is  no  resistance  to  compaction.  The 
solid  Hugoniot  is  vertical  which  means  that  the  end  state  for  compression  is 
always  the  same  volume.  Hence,  the  snowplow  model  is  appropriate  only  to 
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describe  porous  materials  that  completely  crush  when  subjected  to  any  shock 
(see  Figure  6)  [7],  However,  in  reality,  most  porous  materials  offer  some 
resistance  to  compaction,  and  thus  exhibit  some  form  of  partial  compaction. 


As  the  snowplow  model  does  not  allow  for  partial  compaction,  another 
porous  compaction  model  known  as  the  P-a  model  is  more  commonly  used,  as 
to  be  able  to  describe  partially  compacted  states.  W.  Herrmann  initially 
developed  the  irreversible  P-a  compaction  model  in  1969  [8]  to  provide  a 
mechanism  to  describe  the  behavior  of  slightly  distended  material.  Since  then, 
there  has  been  a  successful  extension  of  the  model  to  materials  with  much  larger 
distensions  [9].  A  porosity  parameter,  a,  has  been  introduced  to  allow  for  the 
convenient  separation  of  volume  change  due  to  compression  of  the  material  from 
that  due  to  pore  collapse. 

Porosity,  cc  =  ^ 

where  V  is  the  specific  volume  of  porous  material  and  Vs  is  the  specific  volume  of 
the  corresponding  solid  material  at  the  same  temperature  and  pressure.  Note 
that  porosity,  a,  becomes  unity  when  the  material  is  solid. 
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The  P-a  compaction  model  describes  the  dynamic  compaction  of  a  porous 
material  [9],  as  shown  in  Figure  7.  There  are  several  stages  to  the  compression 
and  release  of  materials  in  this  model; 

a.  Initial  elastic  compression  is  governed  by  elastic  buckling  of  cell 
walls. 

b.  The  beginning  of  permanent  volume  change  coincides  with  the  start 
of  plastic  deformation  of  the  cell  wall  at  pressure  Pe  (which  is  the 
yield  strength  of  the  porous  material) 

c.  Dynamic  compaction  then  moves  along  Rayleigh  lines  to  end  states 
on  the  plastic  curve  (described  by  the  equation)  until  the  pressure 
reaches  the  fully  compacted  state  at  the  solid  compaction  pressure 
Ps  or  is  released  to  a  partially  compacted  state. 

d.  Unloading  from  a  partially  compacted  state  is  elastic  with  its  end 
state  at  a  smaller  ambient  specific  volume  than  the  starting  volume 
due  to  localized  densification. 


Figure  7.  P-a  Compaction  Model 

This  model  is  the  simplest  model  that  can  realistically  describe  the  shock 
compression  and  release  process  in  porous  materials.  This  model  also 
adequately  describes  the  release  from  partially  compacted  states,  which  is 
important  for  low-stress  applications.  Because  of  this,  the  P-a  model  is  widely 
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used  in  the  dynamic  materials  community,  and  is  widely  used  in  hydrodynamic 
computer  code  simulations.  More  advanced  models  exist,  but  they  are  typically 
more  complicated  and  computationally  intensive.  For  this  reason,  the  author  has 
chosen  to  use  the  P-a  model. 
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III.  HYDROCODE  SIMULATIONS 


A.  AUTODYN®  HYDROCODE 

A  hydrocode  is  a  numerical  computer  program  that  uses  a  combination  of 
finite  difference  and  finite  element  techniques  to  solve  dynamic  problems  that 
occur  in  a  short  time  scale  [10],  [11].  Examples  of  such  problems  include  kinetic 
energy  penetration  phenomena,  for  example,  a  projectile  or  bullet  impacting  a 
barrier.  For  this  study,  the  author  uses  the  AUTODYN®2D  hydrocode  (by  Century 
Dynamics)  to  model  the  non-linear  behavior  of  a  projectile  impacting  a  multi¬ 
layered  structure  so  he  can  model  its  ballistic  performance.  He  performs  the 
numerical  computation  in  AUTODYN®  using  the  governing  equations  of 
conservation  of  mass,  momentum  and  energy,  as  defined  in  the  Euler  equations. 
In  addition,  there  must  be  specification  of  the  material’s  equation  of  state  and  a 
suitable  constitutive  strength  model  to  accurately  describe  the  response  of  solid 
materials  and  their  behavior  when  subjected  to  transient  shock  loads.  As  for  all 
numerically  different  computer  codes,  artificial  viscosity  is  used  to  damp  out 
numerical  instabilities  that  can  arise  for  the  particular  differencing  technique 
used.  Although  there  is  a  need  for  artificial  viscosity  to  achieve  stable 
calculations,  this  can  lead  to  slight  errors  in  stress  and  energy. 

B.  GENERAL  HYDROCODE  SETTINGS 

There  are  several  critical  settings  available  in  AUTODYN®  for  setting  up  a 
simulation.  Some  of  the  discussed  settings  will  determine  how  accurately  the 
simulations  model  the  actual  scenario. 

1.  Meshing  Description 

In  hydrocodes,  there  are  two  main  types  of  grid  systems  used  to  describe 
material  movement,  namely  the  Lagrange  and  Euler  system.  The  author  uses 
the  Lagrange  system  in  the  numerical  simulations  in  this  study.  Unlike  an 
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Eulerian  code  which  has  a  fixed  grid,  the  Lagrange  code  has  a  grid  embedded 
with  the  material,  and  as  the  code  tracks  the  flow  or  motion  of  individual  masses, 
the  grid  deforms  together  with  the  material.  Hence,  Lagrange  codes  are  simpler 
and  require  fewer  computations  since  they  do  not  require  a  transport  algorithm  to 
move  material  to  neighboring  cells.  However,  Lagrange  codes  have  problems 
with  accurately  tracking  large  deformations.  With  large  deformations,  numerical 
problems  arising  from  grid  distortion  and  tangling  of  the  mesh  can  lead  to  loss  of 
accuracy  and  time  steps  becoming  small  enough  to  terminate  computation.  To 
overcome  this,  and  especially  for  highly  localized  distortions,  an  erosion 
algorithm  has  to  be  implemented.  The  erosion  algorithm  removes  the  distorted 
grid  elements  when  they  reach  a  predefined  criterion,  which  is  normally  some 
value  of  plastic  strain.  This  allows  the  computation  to  continue,  but  there  is  no 
preservation  of  the  conservation  of  energy  since  there  is  no  longer  tracking  of  the 
internal  energy  of  the  failed  elements.  However,  previous  studies  have  shown 
that  close  agreement  between  calculations  employing  erosion  techniques  and 
experimental  results  can  be  achieved  [12]. 

2.  Simulation  Setup 

The  author  uses  a  2D  axisymmetric  geometry  in  the  calculations  so  that 
he  can  complete  the  simulation  in  a  reasonable  time.  In  numerical  analyses  of 
dynamic  loading,  it  is  necessary  to  use  fine  meshes  to  ensure  the  accuracy  of  the 
results.  However,  simulations  using  short  element  lengths  are  computationally 
intensive  and  take  much  longer  time.  Hence,  the  author  has  to  select  an 
appropriate  mesh  size  that  gives  good  fidelity  but  with  reasonable  run  times.  For 
this  study,  the  author  uses  square  elements  of  0.5  mm  or  two  elements  per  mm 
in  his  simulations. 

In  the  simulations  to  replicate  the  impact  experiments,  the  author  uses 
cylindrical  projectiles  made  of  A2  steel  rods  and  heat-treated  to  a  hardness  of  55- 
58  HRC  (The  hardness  value  of  HRC  55  is  equivalent  to  material  yield  strength 
of  1.8  GPa).  The  projectile  lengths  range  from  %  in  (~19mm)  to  1  in  (~25mm) 
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and  the  diameter  is  0.297  in  (~7.5mm).  The  author  assigns  a  constant  initial 
impact  velocity  to  the  cells  covering  the  boundary  of  the  projectile.  For  target 
plates,  he  sizes  the  models  according  to  the  actual  experimental  samples  used 
with  zero  initial  velocity  constraint  imposed.  He  places  gauges  at  predetermined 
fixed  points  within  the  computation  space  and  along  the  axis  in  both  the  projectile 
and  target  to  measure  dynamic  properties  such  as  pressure  and  material 
velocity.  He  then  uses  the  output  of  these  gauges  to  track  the  shock  wave 
propagation  and  analyze  the  pressure  attenuation  characteristics  in  the  target 
materials. 

3.  Material  Models 

An  equation  of  state  (EOS)  is  a  fundamental  thermodynamic  relationship 
between  pressure,  specific  volume,  and  specific  internal  energy.  Since  the  focus 
of  this  study  is  on  porous  materials,  the  author  will  only  discuss  the  porous  P-a 
equation  of  the  state  compaction  model  option  as  used  in  AUTODYN®  in  detail. 
This  actually  consists  of  both  a  compaction  model  and  an  EOS  for  the  solid 
material,  and  so  it  requires  the  user  to  specify  several  critical  parameters  that  will 
describe  the  crush  behavior  of  the  porous  material.  The  parameters  are  as 
follows: 

a.  Initial  density  of  the  porous  material  po 

b.  Bulk  sound  speed  in  the  elastic  compaction  region 

c.  Maximum  elastic  pressure  (pressure  at  yield)  Pe 

d.  Solid  compaction  pressure  Ps 

e.  Solid  material  Hugoniot  parameters  C  and  S. 

The  initial  density  of  the  porous  material  can  be  determined  experimentally 
using  the  immersion  density  technique  or  more  frequently,  it  is  a  parameter 
provided  as  a  material  specification  by  the  material  supplier. 


17 


The  other  input  parameters  for  the  computation  model  are  determined 
using  the  methodology  discussed  in  Grady  and  Winfree  [13].  First,  to  determine 
the  bulk  sound  speed,  shown  below  is  the  relationship  of  how  bulk  sound  speed 
Ce  varies  with  initial  density  po  (3.1): 

Bulk  sound  speed,  ,  where  a:  is  the  bulk  modulus  (3.1) 

For  isotropic  materials,  the  bulk  modulus  a- is  related  to  the  Young’s 
modulus,  E,  and  Poisson’s  ratio  v  by  (3.2): 

Bulk  modulus,  /c  =  ,  (3.2) 

3(1 -2  v)  '  ' 


Using  the  scaling  relation  (3.3)  developed  by  Gibson  and  Ashby  [3],  [14], 
the  expectation  is  for  the  Young’s  modulus  E  to  vary  with  the  square  of  density. 
The  scaling  relations  for  porous  materials  are  derived  mostly  through  empirical 
fitting  of  experimental  data  from  extensive  testing  on  polymeric  foams. 


E 
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(3.3) 


Assuming  material  isotropy  during  deformation,  the  bulk  wave  speed  of 
the  porous  material  can  be  related  to  the  properties  of  the  fully  dense  solid  by 
(3.4): 


Bulk  sound  speed,  c^ 


EsPo 


'3(1-2v)a 


(3.4) 


In  shock  wave  loading,  the  yield  stress  determines  the  Flugoniot  elastic 
limit  Ohel-  Assuming  linear  elastic  behavior  to  yield  and  a  von  Mises  yield 
condition,  ohel  can  be  determined  from  uniaxial  stress  loading  yield  strength  Oy 
by  the  relation  defined  in  (3.5). 


^HEL  = 


^-2o 


(3.5) 
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Using  the  Qy  data  provided  by  the  foam  supplier,  the  author  calculates  the 
value  of  Ohe-  He  then  determines  the  maximum  elastic  pressure,  Pg  at  which 
yielding  initiates  from  the  equation  below  as  defined  in  (3.6).  The  Poisson  ratio 
o  of  the  material  comes  from  published  literature  and  from  supplier  information 
[14],  [15].  The  hydrostatic  pressure  at  which  yielding  occurs  comes  from  Grady 
and  Winfree  [13]: 


Maximum  elastic  pressure. 


1  +  v 


(3.6) 


where  Ohel  is  known  to  vary  with  foam  initial  density  according  to: 
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(3.7) 


This  model,  developed  by  Gibson  and  Ashby  [14],  predicts  that  Ohel 
increases  with  density  as  defined  by  (3.7). 

The  author  calculates  the  values  of  Ohel  for  different  foam  densities  of  the 
same  material  using  the  Oy  data  provided  by  the  foam  supplier.  The  constant  Cy 
can  then  be  determined  from  a  least-squares  fit. 

Finally,  to  determine  the  pressure  for  complete  compaction,  Ps,  the  author 
uses  the  theoretical  relation  of  Ohel  as  a  function  of  density.  The  solid 
compaction  pressure  is  therefore  the  elastic  pressure  Pg  at  which  Pq=  p^.  He 
obtains  the  material  properties  required  to  compute  the  simulation  inputs  for  the 
P-a  compaction  model  from  material  specifications  from  the  supplier  (for  FR- 
6700  rigid  polyurethane  foam)  and  Ashby  et  al.  [14]  for  the  aluminum  metal  foam. 
Figure  8  shows  the  P-a  compaction  model  of  FR-6720  polyurethane  foam. 

The  solid  material  Hugoniot  parameters,  namely  solid  bulk  sound  speed, 
C  and  the  Hugoniot  shock  velocity  slope  S,  come  from  published  literature  [13]. 
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Figure  8. 


P-a  Compaction  Model  of  FR-6720  Polyurethane  Foam 


In  order  to  completely  describe  porous  material  dynamic  behavior,  the 
author  uses  the  P-a  compaction  model  together  with  the  von  Mises  yield  strength 
criterion  to  model  the  strength  effects.  The  von  Mises  yield  criterion  describes 
the  material  elastic  limit  and  its  transition  to  permanent  plastic  strain  due  to  its 
inability  to  support  large  shear  stresses.  It  is  a  simple  criterion  to  apply,  requiring 
only  material  shear  modulus  and  yield  strength  inputs,  and  is  conveniently 
applied  when  there  is  limited  material  strength  information  available.  Material 
failure  is  modeled  separately,  and  occurs  when  the  material  is  not  able  to 
withstand  tensile  stresses  exceeding  the  material’s  local  tensile  strength.  The 
author  uses  the  hydrodynamic  tensile  (Pmin)  model  for  these  simulations,  and  the 
model  requires  a  specified  constant  hydrodynamic  tensile  limit  to  determine 
failure  occurrence. 

The  author  has  also  chosen  to  use  a  commercially  available  polymeric 
fiber  material,  Dyneema®,  with  very  unique  properties  as  the  wave-spreading 
layer  in  his  armor  concept.  The  Dyneema®  material  is  made  of  several  layers  of 
ultra  high  molecular  weight  polyethylene  (UFIMWPE)  fibers,  with  identical 
perpendicular  fiber  orientation  in  each  layer.  This  orthogonal  configuration  allows 
for  the  distribution  of  the  impact  energy  from  the  penetrator  along  the  fibers  much 
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faster  and  more  efficiently  than  for  an  isotropic  material.  Dyneema®  material  is 
often  used  in  the  design  of  protective  gear  such  as  bulletproof  vests.  The 
orthotropic  material  model  in  AUTODYN®  is  commonly  used  to  describe 
materials  with  different  material  properties  or  strengths  in  different  orthogonal 
directions,  and  in  this  study,  it  will  be  used  to  describe  the  wave-spreading 
characteristics  of  Dyneema®.  Assuming  identical  properties  in  the  fiber 
directions,  the  author’s  chosen  polymeric  fiber  material  falls  in  a  special  class  of 
orthotropic  materials  known  as  “traversely  isotropic.”  As  such,  the  longitudinal 
and  shear  sound  speeds  in  the  two  fiber  directions  in  the  layer  plane  will  be 
identical.  The  sound  speeds  in  the  fiber  directions  will  also  be  much  larger  than 
the  through-thickness  sound  speeds  due  to  its  anisotropy.  This  leads  to  higher 
elastic  wave  speeds  and  Young’s  modulii  E22  and  E33  in  the  fiber  (or  traverse) 
direction  compared  to  Young’s  modulus  En  in  the  through-thickness  (or 
longitudinal)  direction.  The  thesis  work  of  Ong  discusses  the  Dyneema® 
orthotropic  material  model  in  detail  [2]. 

The  other  materials  in  the  armor  concept  use  resident  material  models 
with  modifications  to  known  properties  to  match  the  actual  material  behavior,  (for 
example,  the  Johnson-Holmquist  2  constitutive  model  for  the  ceramics).  Figure  9 
shows  the  material  layering  sequence  for  both  the  numerical  and  experimental 
studies.  Table  1  summarizes  the  material  properties  used  in  the  simulations. 


Projectile 
(Steel  A2) 


□ 


High  Impedance  Layer  (Alumina, 
AI2O3  based  ceramic) 

I  I  Wave-spreading  Layer  (Dyneema®) 

Porous  Layer  (Porous  Foam) 


■  Final  Stopping/Confinement  Layer 
(AI6061-T6) 


Figure  9.  Material  Layering  Sequence  for  Simulations 
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Material 

Equation  of  State 

Constitutive 
Strength  Modei 

Failure  Model 

Steel  A2 

Shock 

C  =  4569  m/s 

S  =  1.49 

Johnson  Cook 

G  =  81.8  GPa 

Ov  =  1.8  GPa 

Hydro  (Pmin) 
r^spaii  — 2  GPa 

Steel  4140 

Linear 

K  =  159  GPa 

Johnson  Cook 

G  =  77  GPa 

Oy  =  0.65  GPa 

Johnson  Cook 
D1=0.05,  D2=3.44 
D3=-2.12,D4=  0.002, 
D5  =  0.61 

Ceramic 

AI203-99.5 

Polynomial 
=  231  GPa 

A2  =  -160  GPa 

A3  =  2774  GPa 

Johnson-Holmquist 

G  =  152  GPa 

Ohel  “  6.57  GPa 

Johnson-Holmquist 
o^spaii  — 0.262  GPa 

Dyneema® 

Din 

Ortho 

Eii  =2.67  GPa 

E22  “  E33  =  96  GPa 
^12  “  ^23  “  0.07 
1/31=0.7 

Gi2  —  G31  =  1  GPa 
G23  =  36  GPa 

Orthotropic  Yield 

A-|  1  =  1 .5 

A22  =  A33  =  1 

Ai3  =  Ai3  =  -0.68 

A23  =  0.26 

A44  =  A55  =  Aee^  4 

Orthotropic  Softening 
aii=  2MPa 

CT22=CT33  =  2000  MPa 

T2  “^23  “^31  MPa 

£11=5.447  KJ/m^ 
£22=^33  =0.3  KJ/m^ 
£12  =  £23=  ^31 

=14.613  KJ/  m^ 

PR-6710 

P-a 

P porous~0.16  g/cm 

Pe  =  2.60  MPa 

Ps  =  112.54  MPa 

C  =  2490  m/s 

S  =  1.56 

von  Mises 

G  =  19.36  MPa 

Oy  =  2.60  MPa 

Hydro  (Pmin) 
r^spaii  — 2  GPa 

Alulight®500 

P-a 

P porous“0.5  g/cm 

Pe  =  7  MPa 

Ps  =  133  MPa 

C  =  5240  m/s 

S  =  1.40 

von  Mises 

G  =  1.88  GPa 

Oy  =  7  MPa 

Hydro  (Pmin) 

Ospaii  =  -2  GPa 

AI6061-T6 

Shock 

Cl  =  5240  m/s 

SI  =  1.4 

Steinburg  Guinan 

G  =  27.6  GPa 

Ov  =  0.29  GPa 

Hydro  (Pmin) 
o^spaii  — 2  GPa 

Table  1.  Material  Properties  for  Numerical  Simulation  in  AUTODYN® 
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c. 


CODE  VALIDATION 


1.  Symmetric  Impact:  Tantalum-Tantalum  Plate  Impact 

The  author  performs  a  simple  plate-to-plate  impact  simulation  to  better 
understand  the  impact  physics  phenomena  and  to  validate  the  AUTODYN® 
hydrocode.  Both  plates  are  of  the  same  material,  namely  tantalum  metal, 
because  Hugoniot  information  for  this  material  is  readily  available  from  the 
literature.  Since  this  a  simple  problem  that  can  be  solved  analytically,  the 
derived  solutions  can  then  be  compared  to  the  simulation  results.  If  the  stress 
wave  propagation  characteristics  are  in  agreement  with  analytical  results,  it  will 
then  give  the  author  confidence  that  the  simulation  setup  is  correct. 

The  2D  symmetric  impact  simulation,  performed  using  the  Lagrange  code 
solver,  consists  of  a  tantalum  flyer  plate  that  is  100  mm  by  5  mm  thick  and  a 
target  plate  of  the  same  material  that  is  100  mm  by  10  mm  thick.  The  author 
uses  a  square  grid  meshing  of  two  elements  per  mm  and  assigns  a  constant 
initial  impact  velocity  of  500  m/s  to  the  flyer  plate.  He  uses  gauges  that  are  set 
up  at  an  interval  distance  of  2  mm  within  the  target  plate  to  measure  impact 
pressure  as  a  function  of  propagation  distance.  Similarly,  he  places  gauges  at  a 
1  mm  interval  distance  within  the  flyer  plate  along  the  axis  of  penetration  (see 
Figure  10). 

The  simulated  shock  pressure  profile  indicated  a  peak  impact  pressure  of 

15.6  GPa  at  the  first  gauge  location  within  the  target  (see  Figure  11).  The 

simulation  results  also  revealed  shock  wave  profiles  that  show  flat  top  waves  with 

short  shock  rise  time  (see  Figure  12).  The  author  observed  a  more  gradual 

release  with  a  decrease  in  pressure,  as  is  expected  for  release  fans.  Using  the 

momentum  jump  condition,  material  properties  (see  Table  2),  and  symmetric 

impact  characteristic  that  the  particle  velocity  Up  in  the  target  is  half  the  flyer  initial 

velocity  Ud,  the  impact  pressure  is  analytically  determined  to  be  15.5  GPa.  The 

small  difference  in  results  and  the  small  wave  perturbations  observed  in  the 

gauge  pressure  plots  are  due  to  numerical  noise  and  artificial  viscosity.  The 
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usage  of  artificial  viscosity  hydrocodes  allows  the  code  to  handle  discontinuities 
such  as  shock  waves  where  there  are  known  issues  with  numerical  instabilities. 
To  model  projectile  penetration  problems  where  shock  waves  are  present,  the 
natural  viscosity  of  most  materials  is  too  small  to  damp  out  numerical  instabilities 
caused  by  the  very  fast  shock  rise.  Hence,  artificial  viscosity  must  be  added  to 
the  pressure,  which  has  the  effect  of  spreading  the  shock  wave  over  several 
mesh  widths,  and  this  transforms  it  from  a  discontinuity  to  a  steep  stress  gradient 
and  makes  computation  more  manageable  [12]. 

Upon  collision  between  impactor  (flyer)  and  target,  a  compression  wave 
propagates  into  both  the  flyer  and  target  with  an  initial  pressure  intensity  of 
P=PoUsUp,  where  po  is  the  material  density.  Us  is  the  shock  velocity  (a 
characteristic  of  the  material),  and  Up  is  the  particle  velocity.  Figures  11  and  12 
show  initial  impact  pressures  of  15.6  GPa  in  both  the  flyer  and  target  plate,  which 
is  consistent  with  the  fundamental  concept  that  if  a  discontinuity  in  pressure  or 
particle  velocity  occurs  across  an  interface,  shock  waves  will  be  created  to  force 
the  pressure  and  Up  to  be  equalized. 

The  analytical  x-t  and  P-Up  diagrams  (see  Figure  13)  also  show  material 
behavior  under  shock  similar  to  the  simulation  results.  For  the  x-t  diagram, 
pressure  and  Up  are  conserved  across  the  impact  surface  at  state  1 .  Since  the 
flyer  plate  is  thin  in  comparison  to  the  target,  the  shock  that  was  created  at 
impact  is  reflected  from  the  backside  free  surface  of  the  flyer  plate  as  a  release 
wave,  resulting  in  a  complete  release  to  zero  pressure  at  state  2.  This  is  similarly 
observed  at  3  when  the  shock  eventually  reaches  the  free  surface  of  the  target 
plate  at  state  3.  Negative  pressure  is  then  created  at  state  4  when  the  two 
release  waves  intersect,  causing  the  target  plate  to  see  localized  tension.  If  the 
dynamic  tensile  strength  of  the  target  material  is  exceeded,  spall  will  occur,  and 
the  target  plate  will  separate  into  two  pieces.  In  Figure  12,  it  can  be  seen  that  at 
a  time  slightly  greater  than  4  ps  the  stress  becomes  negative.  This  corresponds 
to  state  4  in  the  x-t  diagram  as  discussed  above,  where  the  two  release  waves 
interact  and  cause  tension. 
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Figure  10.  Simulation  Setup  for  Tantalum-Tantalum  Plate  Impact 
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Figure  11.  Shock  Pressure  Profile  of  Tantalum-Tantalum  Plate  Impact 


25 


PRESSURE  (kPa ) 


(a)  (b) 

Figure  12.  Pressure  versus  Time  Plot  of  (a)  Target  Plate;  (b)  Flyer  Plate 


Figure  13.  x-t  and  P-Up  Diagram  for  Symmetric  Plate  Impact 
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This  negative  pressure,  or  dynamic  tension,  is  observed  in  the  pressure 
wave  profiles  shown  in  Figure  12(a)  for  the  target  plate  for  several  gauge 
locations. 

To  determine  the  linear  Us-Up  for  tantalum,  the  author  performed 
simulations  using  various  initial  impact  velocities  for  the  impactor.  The  impact 
pressure  is  then  determined  from  the  pressure  plots,  and  by  applying  the 
momentum  jump  condition  the  shock  velocity  Us  is  calculated.  The  linear  Us-Up 
determined  through  simulations  was  found  to  be  Us=1 .38up+3.30  (see  Figure  14), 
which  agrees  reasonably  well  with  the  material  properties  of  tantalum  in  Table  2. 
The  author  calculates  this  as  a  check  to  make  sure  the  computer  code  is 
performing  as  expected. 


Figure  14.  Linear  Us-Up  for  Tantalum  Metal 


Since  the  simulation  results  are  both  qualitatively  and  quantitatively  similar 
to  analytical  solutions,  it  gives  the  author  confidence  that  the  simulation  setup 
accurately  models  the  material  behavior  during  impact. 
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2.  Unsymmetrical  Impact:  Tungsten  Projectile  Impacting 
Tantalum  Plate 

To  validate  the  simulation  setup  for  an  unsymmetrical  projectile  impact, 
the  author  impacted  a  rod-like  tungsten  projectile  (4  mm  in  diameter)  with 
constant  velocity  Ud  of  500  m/s  against  a  large  flat  tantalum  target  plate  that  is 
100  mm  X  100  mm  by  15  mm  thick.  He  used  gauges  that  were  set  up  at  an 
interval  distance  of  2  mm  within  the  target  plate  to  measure  impact  pressure  as  a 
function  of  propagation  distance.  Similarly,  he  placed  gauges  at  a  2  mm  interval 
distance  within  the  projectile  along  the  axis  of  penetration  (see  Figure  15). 

Upon  impact,  compressive  shock  waves  propagated  into  both  the  target 
and  the  projectile,  emanating  from  the  interface.  Figure  16  shows  the  edges  of 
the  planar  waves  in  the  target  material  being  nibbled  away  within  a  short  time 
due  to  lateral  release  waves  from  the  free  surface  where  contact  with  the 
projectile  ceased.  This  happens  in  this  geometry  much  more  quickly  than  for  the 
simulation  above,  which  used  a  flyer  plate  of  the  same  diameter  as  the  target. 
These  lateral  release  waves  reduce  the  intensity  of  the  compressive  waves  to 
form  the  spherical  compressive  pulse,  and  cause  the  initial  state  of  uniaxial  strain 
to  become  that  of  two-dimensional  strain.  During  the  simulation,  the  author 
observed  compression  of  the  ductile  tantalum  target  plate  in  both  the  longitudinal 
and  radial  directions  by  the  penetrating  tungsten  projectile.  Since  the  lead  shock 
is  being  relieved  from  the  sides  as  it  moves  through  the  target,  it  quickly  takes  on 
a  triangular  wave  profile.  This  occurs  when  release  waves  reach  the  shock  front, 
and  indicates  that  the  shock  is  no  longer  supported  from  behind  by  the  projectile. 
When  the  triangular  compressive  wave  eventually  reaches  the  back-free  surface 
of  the  target  plate,  it  is  reflected  from  the  free  surface  as  a  release  (or  tensile) 
wave.  The  negative  pressure  (see  Figure  17)  observed  in  the  pressure  profile  in 
the  tantalum  target  is  the  result  of  this  wave  reflection.  This  causes  created 
tension  within  the  target  material.  If  the  magnitude  and  duration  of  the  tensile 
stress  waves  are  sufficient,  material  failure,  resulting  in  physical  separation  of  the 
material  in  a  variety  of  failure  mechanisms  can  occur.  For  a  ductile  material  like 
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tantalum,  this  failure  is  mainly  due  to  ductile  void  growth,  nucleation,  and 
coalescence.  The  simulation  results  also  confirmed  that  triangular  shock  waves 
are  quickly  obtained  in  the  target,  unlike  the  flat,  top  shock  waves  obtained  for 
flat  plate  impacts,  which  results  in  their  attenuation  as  they  move  through  the 
target.  This  is  expected  as  the  small  projectile  experiences  edge  releases  much 
faster  compared  to  a  large  flat  plate. 
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Figure  15.  Simulation  Setup  for  Tungsten  Projectile  Impacting  Tantalum  Plate 


Figure  1 6.  Pressure  Plot  at  Initial  Impact  after  (a)  0.5  ps;  and  (b)  1 .5  ps 
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TIME  (ms) 

Figure  17.  Pressure  Versus  Time  Plot  of  Tantalum  Target  Plate 

From  the  AUTODYN®  simulation  results,  a  peak  pressure  of  18.2  GPa 
was  obtained.  Using  a  technique  known  as  impedance  matching  for 
unsymmetrical;  impact,  the  peak  pressure  can  be  analytically  determined  by 
using  the  material  properties  from  the  AUTODYN®  material  library  (see  Table  2) 
and  applying  equal  pressure  and  particle  velocity  boundary  conditions  at  the 
interface. 


Material 

Tungsten 

Density,  poi 

19.235  gW 

Linear  Us-Ud  relationship 

U3  =  4.04+1  .23ud  ,  where  Ci=  4.04,  Si=1 .23 

Material 

Tantalum 

Density,  poi 

16.656  g/cm'^ 

Linear  Us-Ud  relationship 

Us  =  3.43+1  .19ud  ,  where  62=  3.43,  82=1 .19 

Table  2.  Material  Properties  from  AUTODYN  Material  Library 
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By  applying  impedance  matching,  Pi  =  P2  (Pi:  pressure  in  projectile;  P2: 
pressure  in  target  plate)  and 

(P01S1-P02S2)  Up2^  -  (2poiSiUD+pOlCi+p02C2)Up2  +  poi(SiUd^  +CiUd)  -  0 
Solving  the  quadratic  equation  derives  Up2  =  0.2844  km/s. 

Then,  P2  =  po2(S2Up2+C2)  Up2^  —  17.85  GPa. 

Checking,  Pi  =  poiSi(Ud-Up2)^  poiCi(Ud"Up2)  —  17.85  GPa. 

From  the  above,  one  can  see  that  the  simulation  results  agree  reasonably 
well  with  the  analytical  results.  Hence,  it  gives  the  author  confidence  that  the 
simulation  setup  can  accurately  model  unsymmetrical  impact  for  a  rod-like 
projectile  impacting  against  a  target  plate. 

3.  Symmetric  Impact:  Porous  Material  (FR-6720  Polyurethane 
Foam) 

Using  the  p-a  material  model  as  the  author  has  defined  it  for  the  FR-6720 
rigid  polyurethane  foam,  he  performed  a  simple  symmetric  impact  simulation  to 
investigate  the  material  properties  and  wave  propagation  characteristics  of  the 
porous  material.  Two  100  mm  square  plates  were  generated  using  planar 
symmetry  in  AUTODYN®.  The  stationary  target  plate  is  20  mm  thick  and  the 
author  assigned  the  10  mm  flyer  plate  an  impact  velocity  of  500  m/s.  Gauges  are 
set  up  at  an  interval  distance  of  2  mm  within  both  the  flyer  and  target  plate  along 
the  axis  of  penetration  to  measure  the  impact  pressure  (see  Figure  18). 

Figure  19  shows  an  initial  impact  pressure  of  0.04  GPa,  which  is  very  low 
compared  to  the  impact  pressure  experienced  by  a  similar  tantalum-tantalum 
impact  simulation.  The  results  also  reveal  that  the  porous  material  has  good 
shock  wave  attenuating  properties,  which  makes  it  suitable  for  energy  absorption 
applications.  However,  due  to  its  low  density  and  hence  poor  resistance  to 
penetration,  it  is  not  suitable  as  a  first  or  second  layer  of  a  multi-layered 
composite  armor  structure.  Hence,  the  author  uses  the  porous  layer  as  the 
energy  absorbing  third  layer  in  his  layered  structure. 
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Figure  18. 


Simulation  Setup  for  Porous  FR-6720  Symmetric  Impact 
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Figure  19. 


Pressure  Versus  Time  Plot  of  FR-6720  Target  Plate 


33 


Detailed  analysis  of  the  shock  profile  in  Figure  19  revealed  multiple 
shocks  propagating  within  the  porous  target  plate.  The  two-wave  structure 
consists  of  an  elastic  wave,  followed  by  a  slower  plastic  wave.  Figure  20 
explains  this  wave  propagation  phenomenon. 


Figure  20.  Elastic-plastic  Transformation  of  Porous  Material.  After  [7] 

If  the  impact  pressure  does  not  exceed  Pi  (the  Flugoniot  elastic  limit  Ohel), 
a  single  elastic  wave  will  propagate  in  the  material.  However,  if  the  magnitude  of 
the  applied  stress  exceeds  Ohel  but  is  not  sufficient  to  reach  state  3,  two  waves 
will  propagate  through  the  medium,  the  elastic  wave  and  the  plastic  wave.  The 
plastic  wave  moves  with  a  speed  that  is  a  function  of  the  slope  of  the  Rayleigh 
line.  It  is  noted  that  the  Ohel  value  of  11.2  MPa  from  the  pressure  plot 
corresponds  reasonably  well  with  analytical  aHEL=15.6  MPa,  which  is  computed 
based  on  Pe  =9.64  MPa  for  the  porous  FR-6720  polyurethane  foam. 

During  the  shock  compression  process  from  the  initial  rest  condition  at 
state  0  to  state  2,  the  shock  wave  first  drives  the  material  to  state  1,  causing  an 
elastic  wave  to  form,  followed  by  a  second  Rayleigh  line  to  state  2,  forming  a 
slower  plastic  wave.  Since  the  gradient  of  the  Rayleigh  line  from  state  1  to  2  is 
less  than  the  Rayleigh  line  gradient  from  state  0  to  1 ,  the  plastic  wave  travels 
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more  slowly  than  the  elastic  wave.  From  the  pressure  wave  profiles  in  Figure  19, 
the  elastic  wave  speed  is  found  to  be  1,100  m/s,  which  agrees  reasonably  well 
with  the  longitudinal  wave  speed  of  1,022  m/s  computed  from  the  longitudinal 
modulus  F  of  the  material.  The  slower  plastic  wave  speed  is  found  to  be  540  m/s 
from  the  pressure  wave  profiles. 

If  the  impact  pressure  is  increased  to  P3,  the  elastic  wave  is  overdriven 
and  a  single  wave  will  propagate  in  the  material.  In  this  wave,  yielding  is  still 
occurring  but  one  cannot  detect  it  from  wave  profiles  because  the  plastic  wave 
speed  has  exceeded  the  elastic  wave  speed.  The  effect  of  increasing  impact 
pressure  is  demonstrated  by  increasing  the  impact  velocity  of  the  flyer.  From  the 
results  shown  in  Figure  21,  one  can  see  that  the  plastic  wave  catches  up  with  the 
elastic  wave  and  will  gradually  disappear  as  impact  velocity  (or  impact  pressure) 
increases.  A  single  wave  will  eventually  form  and  this  phenomenon  can  be 
verified  through  experiments. 
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Figure  21 .  Plastic  Waves  Generated  at  Impact  Velocity  of  1 ,000  m/s 

Since  the  above  results  showed  credible  shock  wave  characteristics  and 
reasonable  shock  pressures,  it  gives  the  author  confidence  that  porous  material 
behavior  has  been  modeled  reasonably  well  using  the  P-a  compaction  model  in 
AUTODYN®. 


35 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


36 


IV.  MATERIAL  CHARACTERIZATION 


Porous  materials  are  characterized  structurally  by  their  cell  topology  (open 
cells,  closed  cells),  relative  density,  cell  size,  and  shape  and  cell  size  distribution 
[14],  This  section  describes  the  characteristics  of  the  two  foam  materials  studied 
in  this  research,  namely  the  aluminum  metal  foam  and  rigid  polyurethane  foam, 
which  the  author  used  in  two  different  initial  densities. 

A.  ALUMINUM  METAL  FOAM 

The  aluminum  porous  material,  supplied  by  Alulight®  of  America  L.P., 
comes  in  the  form  of  an  aluminum  foam  sandwich.  The  aluminum  foam 
sandwich  (AFS)  is  a  flat  structure  composed  of  a  closed-cell,  foamed  aluminum 
core,  and  aluminum  cover  sheets.  The  sandwich  panel  is  formed  by  foaming  a 
powder  metallurgy-precursor  between  the  metal  cover  sheets.  The  powder 
metallurgy  process  involves  blending  atomized  aluminum  alloy  powder  with 
titanium  hydride  (TiH2)  and  compacting  the  mixture  to  form  a  precursor  material. 
The  compacted  precursor  is  then  placed  inside  a  mold  and  heated  to  a 
temperature  near  the  melting  point  of  the  aluminum  alloy.  The  foaming  agent 
decomposes  when  subjected  to  heat  and  gas  bubbles  that  get  trapped  form 
within  the  compacted  powder  body.  The  expanding  gas  bubbles  then  create 
voids  within  the  expanding  molten  metal  and  the  mold  is  rapidly  cooled  to  prevent 
collapse  of  the  foam  structure.  During  the  expansion  process,  the  liquid  foam  is 
also  metallurgically  bonded  to  the  solid  cover  sheets.  Upon  solidification,  the 
cellular  metal  foam  structure  is  formed  [16]. 

Prior  to  physical  examination  of  the  foam  structure,  specimens  of  size  20 

mm  X  20  mm  were  cut  from  the  large  panel  of  500  mm  x  500  mm  of  thickness  12 

mm.  The  AFS  consists  of  a  10  mm-thick  foam  metal  core  sandwiched  by  1  mm 

thick  aluminum  cover  sheets.  To  avoid  damage  to  the  examined  foam  surface,  a 

slow-speed  diamond  cutter  was  used  to  cut  the  specimen  to  size.  Figure  22 

shows  the  close-up  view  of  a  typical  specimen.  Using  an  open-source  image 
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analysis  software,  Imaged,  parameters  such  as  mean  pore  area,  equivalent  pore 
diameter,  area  fraction,  and  pore  size  distribution  were  determined.  Table  3 
shows  a  summary  of  the  material  characteristics  and  Figure  23  shows  the 
graphical  representation  of  pore  size  distribution  of  the  metal  foam. 


Material 

Alulight®500 

Foam  type 

Closed  cell 

Nominal  density  of  foamed 
metal  core  (g/cm^) 

0.5 

Relative  density 

0.185 

Mean  pore  size  (mm) 

2.27 

Maximum  pore  size  (mm) 

4.56 

Minimum  pore  size  (mm) 

0.24 

Estimated  pore  volume  (%) 

75 

Table  3.  Summary  of  Material  Characteristics  of  Alulight®500  Metal  Foam 


Visual  inspection  of  the  metal  foam  sections  revealed  that  the  foam  has 
non-uniform  pore  structure  with  variable  pore  size.  These  effects  can  be 
explained  by  the  unequal  heating  rate  during  manufacturing  which  leads  to 
variable  pore  size.  Each  cell/pore  is  separated  from  another  cell  by  a  thin  wall  of 
metal.  The  results  showed  that  the  pore  diameter  ranges  from  0.24  to  4.56  mm 
with  the  majority  of  pores  within  a  size  range  of  1  to  4  mm  in  diameter.  The  size 
and  distribution  of  pores  within  metal  foams  are  usually  random.  The  computed 
pore  volume  fraction  of  75%  is  also  representative  of  a  typical  porosity  level  of 
70-90%  for  powder-metallurgy  metal  foams  [17]. 
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Peicentage 


Figure  22.  Close-up  View  of  Typical  Alulight®  AFS  Specimen 


Pofg  Si?g  Distibution  Q^  AluligfctSOO 


Figure  23.  Pore  Size  Distribution  of  ALULIGFIT®500  Metal  Foam 
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B.  RIGID  POLYURETHANE  FOAM 


The  rigid  polyurethane  foam,  supplied  by  General  Plastics  Mfg  Co.,  comes 
in  sheet  form  5  mm  in  thickness.  Rigid  polyurethane  (PU)  foams  are  commonly 
used  in  composite  constructions,  particularly  as  honeycomb  aircraft-interior 
panels  because  of  their  low  weight  and  high  temperature  compressive  strength 
properties.  Polyurethane  foams  are  usually  manufactured  by  a  continuous- 
extrusion  process  to  form  blocks  that  are  then  cut  to  make  sheets.  The  FR-6700 
rigid  polyurethane  foam  is  a  closed-cell,  flame  retardant  polyurethane  foam 
available  in  densities  ranging  from  10  to  25  Ibs/ft^  (0.16  to  0.4  g/cm^).  In  this 
study,  the  author  only  investigates  two  foam  densities;  FR-6710,  10  Ibs/ft^  (0.16 
g/cm^)  and  FR-6720,  20  Ibs/ft^  (0.32  g/cm^). 

Specimens  of  size  50  mm  x  50  mm  were  cut  and  examined  under  an 
optical  microscope.  Figure  24  and  Figure  26  show  microscopic  images  of  the 
FR-6710  and  FR-6720  PU  foam  respectively.  Using  the  same  image  analysis 
software,  the  author  determined  the  microstructural  parameters.  Figure  25  and 
27  show  the  graphical  representation  of  pore  size  distribution  of  the  PU  foam. 
Table  4  summarizes  the  material  characteristics  of  the  PU  foam. 


Material 

FR-6710 

FR-6720 

Foam  type 

Closed  cell 

Closed  cell 

Nominal  density  (g/cm'^) 

0.16 

0.32 

Relative  density 

0.126 

0.253 

Mean  pore  size  (pm) 

133.6 

133.4 

Maximum  pore  size  (pm) 

266.9 

243.8 

Minimum  pore  size  (pm) 

63.4 

59.8 

Estimated  pore  volume  (%) 
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Table  4.  Summary  of  Material  Characteristics  of  Polyurethane  Foams 


Microscopic  examination  of  the  polyurethane  foams  revealed  that  the 

foam  has  variable  spherical  pores  with  diameters  ranging  from  50  to  300  pm. 

Both  foam  densities  exhibit  similar  characteristics  of  mean  pore  diameter  and 

pore  size  distribution,  indicating  that  the  foam  manufacturing  process  is  likely  to 

be  identical.  The  author  also  observed  that  the  pore  volume  of  the  lower  density 
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foam  FR-6710  is  almost  twice  of  that  of  the  higher  density  foam.  This  correlates 
well  with  the  nominal  foam  density  values  claimed  by  the  supplier. 
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Figure  24.  Microscopic  Image  of  FR-6710  Polyurethane  Foam 


Pot»  Si2^  Dislibution  oi  FR-6710  Rigid  Polfut^than^  Foam 


Figure  25.  Pore  Size  Distribution  of  FR-6710  Polyurethane  Foam 
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Figure  26.  Microscopic  Image  of  FR-6720  Polyurethane  Foam 


Polf  Sizt  Distibution  ol  FR-6720  Rigid  Polpuigthang  Foam 


Figure  27.  Pore  Size  Distribution  of  FR-6720  Polyurethane  Foam 
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V.  EXPERIMENT 


A.  INTRODUCTION 

Since  the  focus  of  this  research  is  to  understand  the  effect  of  having  a 
layer  of  porous  material  in  an  advanced  layered  armor  concept,  it  is  appropriate 
to  do  integral  tests  as  a  final  validation.  Since  this  layer  is  one  of  several,  the 
author  developed  a  test  matrix  for  needed  integral  penetration  experiments  in 
close  coordination  with  Ong  [2],  who  has  focused  on  the  orthotropic  material 
layer.  Since  the  author  was  not  able  to  measure  the  fundamental  dynamic 
properties  of  the  foam  materials,  he  used  literature  values  as  input  into  the 
chosen  material  models  used  in  the  computer  simulations.  This  allowed  them  to 
calculate  the  system  response  of  each  of  the  layered  target  geometries,  and 
Mowed  for  a  close  comparison  with  experimental  results.  In  order  to  test  the 
effectiveness  of  the  porous  materials,  the  author  included  what  he  believed  to  be 
appropriate  tests  with  and  without  this  layer.  This  section  documents  the 
experimental  results. 

B.  TEST  SAMPLE  PREPARATION 

The  materials  used  to  fabricate  the  layered  structure  were  cut  to  the 
square  test  sample  size  of  100  mm  by  100  mm.  The  mating  surfaces  were  then 
cleaned  thoroughly  with  alcohol  to  remove  surface  contamination.  The  materials 
were  then  bonded  together  using  a  low  viscosity  epoxy,  AngstomBond  9110LV. 
This  two-part  adhesive  exhibits  good  wetting,  cures  at  room  temperature,  and 
develops  strong,  low-shrinkage  bonds  to  most  materials  including  glass 
ceramics,  metal,  and  rigid  plastics.  The  author  selected  a  slow-curing  and  low- 
viscosity  adhesive  to  minimize  trapped  air  bubbles  between  the  bond  layers, 
which  may  affect  the  ballistic  performance  of  the  composite  layered  structure. 
After  sufficient  application  of  the  adhesive  to  the  mating  surfaces,  the  author 
applied  a  controlled  clamping  force  to  the  bonded  layers  to  purge  the  excess 
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adhesive  so  that  a  thin  bond  layer  could  be  formed  (see  Figure  28).  The 
adhesive  was  then  allowed  to  cure  for  24  hours  under  a  constant  applied  load. 
Extra  care  was  taken  to  ensure  that  the  applied  load  did  not  crush  the  pores  of 
the  porous  materials. 


(a)  (b) 

Figure  28.  Test  Sample  Preparation 

C.  EXPERIMENTAL  SETUP 

Impact  experiments  were  conducted  on  a  compressed-gas  gun  at  the 
facilities  of  University  of  California,  Santa  Barbara  (UCSB).  With  a  maximum  gas 
pressure  of  2,000  psi,  the  maximum  attainable  projectile  velocity  for  a  0.3” 
diameter  steel  rod  of  1”  length  was  found  to  be  484  m/s.  Fligher  projectile 
velocities  can  be  achieved  by  using  smaller  and  lighter  projectiles.  Photographs 
of  the  experimental  set-up  are  given  in  Figure  29. 
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Figure  29.  Photograph  of  Compressed-Gas  Gun  Facility 


The  gas  reservoir  and  the  barrel  were  connected  via  a  flange  assembly. 
The  reservoir  and  the  barrel  could  be  separated  by  removing  the  securing  bolts, 
and  this  enabled  the  projectile  to  be  loaded  into  the  7.62  mm  caliber  smooth  bore 
barrel.  A  steel  diaphragm  was  then  inserted  between  the  flange  interfaces 
before  the  flange  assembly  was  re-assembled.  When  the  reservoir  was 
pressurized  up  to  about  1,800  psi,  the  helium  gas  pressure  ruptured  the  steel 
diaphragm  and  propelled  the  projectile  towards  the  target,  which  was  held  in  the 
impact  chamber.  The  impact  event  was  captured  using  a  high-speed  camera 
viewed  through  the  transparent  armor-glass  in  the  impact  chamber.  The  high- 
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speed  camera  used  in  these  experiments  to  capture  images  at  about  50,000  fps 
was  triggered  by  a  trip-wire  detector  that  was  placed  at  the  exit  of  the  barrel 
muzzle.  Delays  for  the  flash  were  calculated  based  on  the  expected  projectile 
velocity.  The  projectile  impact  velocity  was  estimated  using  image  analysis 
software  to  analyze  the  sequence  of  still  picture  frames  captured  by  the  high¬ 
speed  camera.  A  projectile  capture  system  using  lightly  compacted  sand  was 
placed  behind  the  target  in  the  event  that  the  projectile  penetrated  the  target. 

D.  EXPERIMENTAL  RESULTS 

Sixteen  high-velocity  impact  tests  with  various  experimental  configurations 
were  conducted  at  the  gas-gun  facility.  As  the  focus  of  this  study  is  on  porous 
materials,  only  experimental  data  beneficial  to  this  study  are  shown  in  Table  5. 
Other  experimental  data  may  be  obtained  from  Ong  [2].  Three  types  of  porous 
materials  were  used  in  this  investigation:  rigid  polyurethane  foam  FR-6710  and 
FR-6720  and  Alulight®500  aluminum  metal  foam  (in  the  form  of  AFS). 
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Shot 

No. 

Test 

Sample 

Areal 

Density 

(g/cm^) 

Projectile 

Description 

Impact 

Velocity 

(m/s) 

Complete 

Penetration 

(Y/N) 

6 

S 

3.76 

Steel  Cylindrical  Rod 

500 

N 

9 

S-P2-AI 

4.32 

Steel  Cylindrical  Rod 

521 

N 

10 

S-P2 

3.92 

Steel  Cylindrical  Rod 

509 

N 

11 

S 

3.76 

1”  Steel  Cylindrical  Rod 

484 

Y 

12 

C-D-M 

3.87 

1”  Steel  Cylindrical  Rod 

465 

N 

13 

C-D-P2-AI 

3.38 

1”  Steel  Cylindrical  Rod 

465 

N 

14 

C-D-P1-AI 

3.30 

1”  Steel  Cylindrical  Rod 

481 

N 

18 

C-D 

2.82 

1”  Steel  Cylindrical  Rod 

483 

N 

Legend: 


Symbol 

Material 

Thickness 

(mm) 

Al 

Al  6061-T6  Aluminum 

1.5 

C 

Al  203-99.5  Ceramic 

6 

D 

Dyneema 

5 

M 

Alulight  aluminum  foam  sandwich 
(AFS) 

12 

PI 

FR-6710  rigid  polyurethane  foam 

5 

P2 

FR-6720  rigid  polyurethane  foam 

5 

S 

AISI4140  steel 

4.76 

Table  5.  Summary  of  Experiment  Results 


The  author  performed  a  visual  inspection  of  the  post-impact  samples  and 

measured  critical  dimensions  of  the  crater  profile.  Table  6  contains  a  summary 

of  the  post-impact  test  measurements  on  test  samples.  Figure  30  shows  the 

schematic  of  a  post-impact  target  plate  with  definitions  of  the  measured 

parameters.  During  the  ballistic  testing,  all  laminates  displayed  attributes  typical 

of  armor  under  ballistic  attack  and  will  be  discussed  in  detail  in  Chapter  VI. 

However,  it  should  be  noted  that  the  recovered  steel  projectiles  showed  brittle 

fracture  failure  upon  impact.  This  is  atypical  of  small-arms  ball  ammunition, 

which  has  a  relatively  soft  steel  core  mass  surrounded  by  a  gliding  jacket.  The 

soft  penetrating  mass  is  made  from  280HV  soft  steel,  which  enables  it  to  deform 

when  penetrating  the  target.  Projectiles  used  for  these  tests  shattered  because 

they  were  made  of  hardened  steel,  and  were  therefore  in  a  brittle  condition. 

Therefore,  the  projectile  that  the  author  and  Ong  [2]  used  for  these  tests 
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exhibited  minimal  deformation  and  failed  through  brittle  cracking  rather  than 
plastic  deformation. 


Figure  30.  Schematic  of  Crater  Profile  Indicating  Definitions  of  Measured 

Parameters 


Shot 

No. 

Test 

Sample 

Bulge 

Width 

(mm) 

Overall 

Depth 

(mm) 

Max. 

Crater 

Depth 

(mm) 

Max.  Crater 
Base 
Diameter 
(mm) 

Max. 

Crater 

Diameter 

(mm) 

a 

b 

c 

D 

e 

6 

S 

19.88 

12.01 

2.26 

13.83 

9.15 

9 

S-P2-AI 

72.94 

15.81 

4.94 

13.54 

9.61 

10 

S-P2 

23.86 

11.63 

1.87 

11.93 

8.05 

11 

S 

14.74 

9.25 

9.25 

12.47 

8.44 

12 

C-D-M 

73.61 

33.57 

17.41 

24.04 

24.04 

13 

C-D-P2-AI 

77.51 

32.64 

17.91 

21.58 

21.58 

14 

C-D-P1-AI 

75.16 

35.12 

21.27 

31.76 

31.76 

18 

C-D 

100 

45.97 

39.36 

25.51 

25.51 

*  All  measurements  made  using  digital  vernier  caliper 


Table  6.  Summary  of  Post-impact  Test  Measurements  on  Test  Samples 
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VI.  DISCUSSION  OF  RESULTS 


The  author  modeled  the  impact  experiments  using  the  numerical 
simulation  setup  discussed  in  Chapter  III  to  match  the  actual  measured 
experimental  projectile  impact  velocity.  He  then  compared  the  results  of  the 
analyses  with  actual  experimental  data.  Upon  validation  of  the  models,  the 
author  will  use  AUTODYN®  simulations  to  predict  the  performance  of  the  various 
proposed  multi-layered  composite  armor  structures  normalized  to  the  same 
projectile  velocity.  The  author  performed  the  analyses  to  achieve  the  objectives 
stated  below. 

A.  VALIDATE  SIMULATION  SETUP  THROUGH  EXPERIMENT 

Figure  31  shows  the  results  of  shot  number  6,  which  was  a  steel  rod 
projectile  impacting  an  AISI  4140  steel  target  plate  at  500  m/s.  The  author 
placed  pressure  gauges  at  1  mm  and  4mm  interval  distances  along  the  axis  of 
penetration  within  the  target  plate  and  projectile  respectively.  When  the  projectile 
and  target  collided,  very  high  pressures  were  generated  in  both  the  projectile  and 
the  target.  For  example,  an  impact  pressure  of  about  9.5  GPa  was  attained  in 
the  target  at  the  location  of  gauge  1 .  The  wave  profiles  for  gauges  1  and  2  have 
relatively  flat  tops,  but  subsequent  gauges  show  triangular  wave  profiles.  The 
triangular  wave  shapes  are  caused  by  release  waves  that  are  formed  at  the  free 
surfaces.  After  2ps,  negative  pressures  are  formed,  indicating  the  presence  of 
intersecting  release  waves.  The  intersecting  release  waves  cause  created 
tension  within  the  target  material,  and  if  both  magnitude  and  duration  are 
sufficient,  material  failure  may  occur.  Approximately  0.1  ms  after  impact,  the 
projectile  was  stopped  with  a  predicted  overall  penetration  depth  b  of  12.5  mm. 
The  localized  plastic  material  deformation  at  the  impact  point  observed  in 
simulation  predictions  was  similar  to  experimental  observations.  As  compared 
with  the  experimental  measurements  for  shot  number  6  in  Table  6  of  the 
experimental  results  section,  there  is  good  agreement  between  the  simulation 
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ABS.VEL  (m/s) 


and  experimental  result.  Simulation  predicted  that  complete  penetration  would 
not  occur  and  the  author  observed  this  experimentally.  Figure  32  shows  the 
post-impact  condition  of  the  target  plate  and  high-speed  photographs  of  the  final 
moments  of  the  impact  sequence  taken  at  30ps  inter-frame  time. 


Target 


Projectile  Velocity  vs 


Figure  31 .  Simulation  Results  of  Shot  Number  6 
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target  plate  arrested 


Figure  32.  Experimental  Results  of  Shot  Number  6 
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Figure  33  shows  the  simulation  results  for  Shot  Number  11,  which  was  a 
1”  Steel  Rod  projectile  impacting  an  AISI  4140  steel  target  plate  at  484  m/s. 
Figure  33  shows  similar  pressure  wave  profiles  as  in  Figure  31,  except  with 
higher  pressures.  Figure  33  also  shows  that  the  1”  projectile  actually  perforates 
the  steel  plate  with  a  residual  velocity  of  100  m/s  at  170  ps  after  impact.  Based 
on  the  velocity-time  plot  of  the  projectile,  target  perforation  occurred  at 
approximately  100  ps  after  impact.  The  author  observed  the  occurrence  of  ductile 
hole  growth  with  significant  plastic  deformation  of  the  target  material.  The 
material  is  pushed  out  of  the  way  of  the  projectile,  revealing  a  through  hole 
approximately  the  size  of  the  diameter  of  the  projectile.  The  simulation  results 
agree  well  with  the  experiment  for  general  observations  and  the  final  crater  size. 
The  high-speed  photography  showed  brittle  failure  of  the  impacting  projectile, 
leading  to  ejected  debris  with  an  estimated  residual  velocity  of  118  m/s,  which  is 
in  good  agreement  with  the  simulation  results.  The  ballistic  test  results  of  Shot 
11  also  showed  that  the  estimated  kinetic  energy  required  for  penetration  of  the 
AISI  4140  steel  plate  is  1,000  joules.  The  author  will  use  the  ballistic 
performance  of  the  AISI  4140  steel  plate  as  a  benchmark  against  the  various 
multi-layered  structure  armor  concepts  he  discusses  in  a  later  part  of  this 
chapter.  Figure  34  shows  the  post-impact  condition  of  the  target  plate  and  high¬ 
speed  photographs  of  the  final  moments  of  the  impact  sequence  taken  at  30  ps 
inter-frame  time. 

From  the  above  simulation  results,  the  author  concludes  the  validation  of 
the  simulation  setup,  projectile,  and  target  steel  material  properties  via  the 
experiment  results. 
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Target 


Projectile  Velocity  vs  Time 
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^  (6)Gauga#  1 1 

(7)Gaugft#  12 
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aw 


Figure  33.  Simulation  Setup  of  Shot  Number  1 1 
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Projectile 


Projectile  impacting  target 


Projectile  debris 


Figure  34.  Experimental  Results  of  Shot  Number  1 1 
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B.  ROLE  OF  INERTIAL  BACKING  PLATE 


Figure  35  shows  the  simulation  results  of  Shot  Number  10,  which  was  a 
3/4”  steel  rod  projectile  impacting  a  S-P2  target  plate  at  509  m/s.  In  this 
experiment,  the  foam  material  has  no  backinp  material,  and  thus  is  unsupported. 
Figure  35  clearly  shows  that  the  compression  shock  wave  is  rapidly  attenuated 
from  4.2  GPa  at  gauge  5  to  23  MPa  at  gauge  8  due  to  the  shock  attenuating 
properties  of  the  porous  material.  This  is  equivalent  to  a  pressure  reduction  of 
99.5%  across  2  mm  of  porous  material.  The  pressure  profiles  in  the  porous  layer 
also  show  the  presence  of  a  two-wave  structure  as  explained  in  Chapter  III. 

Figure  36  shows  the  x-t  diagram  of  the  wave  interactions.  When  the  first 
shock  reaches  the  foam  material  interface,  the  shock  reflects  off  the  interface  as 
a  release  wave  into  the  steel  material  and  transmits  into  the  porous  material  as  a 
shock  travelling  at  a  lower  shock  speed.  The  shock  speed  in  the  porous  material 
is  lower  than  the  shock  speed  in  the  steel  plate  because  the  foam  material  has  a 
lower  shock  impedance  and  because  this  shock  causes  pore  collapse.  Since  by 
this  time  the  shock  wave  profile  is  triangular  (unsupported),  it  attenuates  as  a 
function  of  run  distance  as  the  following  release  wave  continuously  erodes  the 
shock  front.  The  author  also  observed  this  effect  in  the  steel  material.  When  the 
transmitted  shock  arrives  at  the  rear  free  surface,  it  reflects  as  a  release  wave, 
which  further  attenuates  the  peak  shock  pressure  within  the  porous  material. 
This  explains  the  pressure  wave  attenuation  recorded  by  the  gauges  in  the 
simulation. 

The  simulation  results  also  showed  signs  of  spallation,  which  is  consistent 
with  observed  experiment  results.  However,  there  was  no  complete  penetration. 
Figure  37  shows  the  experimental  results.  Spallation  is  a  process  in  which 
fragments  of  material  are  ejected  from  a  body  due  to  wave  interactions  that 
cause  dynamic  tensile  state  in  the  material  [12].  Spall  occurs  when  the  dynamic 
tensile  strength  of  the  material  is  exceeded.  Since  there  is  a  triangular  wave 
shape,  tension  will  be  generated  immediately  upon  reflection  from  the  rear  free 
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surface.  The  author  does  not  expect  this  polyurethane  foam  material  to  have  a 
very  high  dynamic  tensile  strength,  so  it  will  spall  very  easily  as  long  as  it  sees 
tension.  The  author  sees  clear  evidence  for  this  process  in  the  simulation, 
where  material  erosion  is  observed  because  of  the  spallation  process. 
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Figure  35.  Simulation  Setup  for  Shot  Number  10 
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Figure  36.  x-t  Diagram  for  S-P2  Target  Plate 
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Figure  37.  Experimental  Results  of  Shot  Number  10 
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Figure  38  shows  the  simulation  results  of  Shot  Number  9,  which  is  a  3/4” 
steel  rod  projectile  impacting  a  S-P2-AI  target  plate  at  521  m/s.  In  this 
experiment,  the  author  added  a  thin  aluminum  layer  to  the  rear  surface  of  the 
porous  layer  to  provide  support.  Figure  39  shows  that  the  porous  material 
attenuates  the  compression  pressure  from  4.5  GPa  at  gauge  5  to  60  MPa  at 
gauge  8.  The  pressure  reduction  is  about  98.6%  across  2  mm  of  the  porous 
material.  The  x-t  diagram  in  Figure  39  explains  the  higher  peak  pressures 
recorded  by  the  gauges.  The  shock  propagating  in  the  porous  material  is 
reflected  off  the  interface  of  the  higher-impedance  aluminum  inertia  backing  plate 
as  a  re-shock,  which  then  increases  the  gauge  pressure.  The  simulation  results 
show  that  the  aluminum  backing  plate  is  plastically  deformed  due  to  the  impact, 
but  the  porous  rigid  polyurethane  foam  is  kept  relatively  intact.  Complete 
penetration  did  not  occur.  The  above  results  show  good  agreement  between  the 
code  and  the  experiment.  Figure  40  shows  the  post-impact  condition  of  the 
target  plate. 

Comparing  this  result  to  that  above  for  foam  with  no  backing  material,  the 
author  concludes  that  the  inertial  backing  plate  was  necessary  to  constrain  the 
porous  material  and  prevent  it  from  spalling.  This  then  enables  the  porous  media 
to  perform  its  energy  absorption  role  more  effectively  later. 
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Figure  38.  Simulation  Setup  for  Shot  Number  9 
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Figure  39.  x-t  Diagram  for  S-P2-AI  Target  Plate 
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Figure  40.  Experimental  Results  of  Shot  Number  9 
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c. 


MULTI-LAYERED  ARMOR  PLATE  PERFORMANCE 


There  are  several  interdependent  requirements  such  as  cost, 
manufacturability,  and  structural  and  ballistic  performance  in  the  design  of  armor 
structures.  While  the  overall  design  objective  is  to  balance  all  these 
requirements,  one  usually  uses  a  critical  ballistic  performance  parameter  such  as 
dynamic  deflection  to  control  the  design.  Dynamic  deflection  is  the  maximum 
deflection  the  armor  undergoes  during  a  ballistic  impact  event.  This  distance  is 
an  indication  of  the  standoff  distance  required  to  avoid  physical  trauma  to  the 
soldier  wearing  the  personnel  armor.  For  vehicle  armor,  the  standoff  distance 
affects  the  space  available  for  loads  and  soldiers  inside  the  vehicle.  On  one 
hand,  large  dynamic  deflections  are  desirable  since  they  generally  result  in 
improved  ballistic  efficiency,  but  they  can  lead  to  undesirable  space  requirements 
in  standoff  distance  [18].  Since  the  focus  of  this  current  study  is  on  the  ballistic 
performance  of  the  multi-layered  armor  concept,  the  author  will  use  dynamic 
deflection  as  the  performance  metric. 

As  it  was  not  possible  to  attain  consistent  projectile  impact  velocity  for  all 
impact  tests,  in  order  to  conduct  performance  comparisons  between  different 
multi-layered  armor  structures  the  author  has  to  first  validate  the  material  model 
parameters  determined  in  Chapter  III  by  performing  AUTODYN®  simulations  to 
replicate  the  live  firing  experiments.  Upon  validation  of  the  models,  the  author 
can  then  normalize  the  initial  projectile  impact  velocities  for  the  simulation  models 
to  the  same  impact  velocity  to  allow  a  comparison  study. 

1.  Ceramic-Dyneema®  (C-D)  Target  Plate 

The  target  plate  consists  of  a  6  mm  ceramic  tile  layered  over  a  5vmm 
Dyneema®  material.  Based  on  the  recorded  experimental  data  from  Shot  18, 
the  1”  length  cylindrical  A2  steel  projectile  impacted  the  composite  plate  at  483 
m/s.  Figure  41  shows  the  simulation  results.  The  pressure  profiles  obtained 
from  the  impact  simulation  showed  that  the  peak  pressure  is  rapidly  attenuated 
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by  the  wave-spreading  properties  of  the  Dyneema®  material.  The  pressure  plots 
also  showed  that  the  wave  speed  in  the  through-thickness  direction  is  much 
smaller  than  the  lateral  wave  speed,  as  expected  for  this  material.  The 
composite  plate  resisted  complete  penetration  and  the  detainment  of  the 
projectile  was  about  0.32  ms  after  impact. 

The  role  of  the  ceramic  layer  in  an  impact  event  is  to  deform  the  projectile 
and  strongly  dissipate  its  energy  before  it  reaches  the  Dyneema®  backing  plate. 
This  happens  by  causing  extensive  plastic  deformation  in  the  projectile, 
effectively  turning  kinetic  energy  into  internal  (heat)  energy.  Ceramic  materials 
possess  very  high  compressive  strength  but  are  very  weak  in  tension.  Hence, 
ceramics  are  extremely  brittle  and  susceptible  to  fracture.  Figure  42  shows  the 
post-impact  condition  of  the  target  plate.  The  ceramic  layer  showed  heavy 
cracking  because  of  the  impact.  Radial  cracks  propagated  away  from  the  impact 
point  to  the  edge  of  the  plate,  and  a  large  quantity  of  ceramic  material  had 
pulverized  and  ejected  during  the  impact.  The  large  deformation  in  the 
Dyneema®  material  resulted  in  layer  delamination  within  the  material  and  fiber 
stretching,  both  of  which  cause  kinetic  energy  to  be  absorbed.  The  failure 
mechanism  observed  in  the  simulation  agrees  well  with  the  experimental  results. 
Furthermore,  the  overall  deformation  measurements  are  also  comparable  (see 
Table  7).  Since  the  simulation  results  are  both  qualitatively  and  quantitatively 
similar  to  the  experimental  results,  it  gives  the  author  confidence  that  there  is 
accurate  modeling  of  the  material  behaviors  in  AUTODYN®. 
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Figure  41 .  Simulation  Results  of  Shot  Number  18 
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Figure  42.  Experimental  Results  of  Shot  Number  18 
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2. 


Ceramic-Dyneema®-FR6720-AI  (C-D-P2-AI)  Target  Plate 


Next,  the  author  looks  at  a  ‘complete’  composite  armor  plate.  This  plate 
has  all  the  components  required  to  defeat  the  penetration  process  discussed 
earlier.  The  target  plate  consists  of  a  6  mm  ceramic  layer,  a  5mm  Dyneema® 
layer,  and  a  5  mm  FR-6720  polyurethane  (PU)  foam  0.32  g/cm^  supported  by  a 
thin  Aluminum  6061 -T6  inertial  backing  plate  to  constrain  the  PU  foam.  The 
previous  section  showed  that  the  inertial  backing  plate  was  necessary  to  confine 
the  porous  material  so  that  it  could  perform  its  energy  absorption  role  efficiently. 
As  before,  the  author  used  a  1”  length  cylindrical  A2  steel  projectile.  The 
projectile  impacted  the  target  plate  at  velocity  of  463  m/s,  according  to 
experimental  data. 

The  stress  wave  propagation  across  a  structure  consisting  of  dissimilar 
materials  can  be  extremely  complex  and  may  not  be  easily  understood  using 
fundamental  shock  physics.  The  impact  problem  quickly  develops  into 
complicated  2D  interfacial  wave  interactions  between  dissimilar  materials.  Edge 
releases  of  the  projectile  at  the  impact  interface  further  complicate  this,  which 
destroys  the  state  of  uniaxial  strain.  Therefore,  computer  codes  such  as 
hydrocodes  are  used  to  simulate  complicated  non-linear  dynamic  problems  such 
as  ballistic  impacts  that  cannot  be  solved  analytically.  Hence,  the  performance  of 
the  porous  foam  can  only  be  quantitatively  analyzed  through  its  pressure 
attenuation  properties.  The  simulated  pressure  plots  in  Figure  43  clearly  show 
that  the  compression  shock  wave  was  strongly  attenuated  by  the  porous 
material.  The  pressure  was  reduced  from  0.2  GPa  at  the  last  gauge  point, 
namely  gauge  11  in  the  Dyneema®  layer  to  82  MPa  at  gauge  16  in  the  PU  foam. 
This  is  equivalent  to  a  pressure  reduction  of  59%  across  5  mm  of  porous  foam. 
Compared  to  the  C-D  composite  plate  configuration,  the  incorporation  of  the 
porous  PU  foam  also  increases  the  stress  wave  rise  time  of  the  propagating 
stress  wave  and  significantly  delays  the  stress  wave  propagation.  This 
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effectively  spreads  out  in  time  the  delivery  of  momentum.  The  hydrocode 
simulation  predictions  capture  both  the  widening  rise  time  and  stress  wave 
attenuation  characteristics  of  the  porous  media. 

The  author  also  observed  that  the  plastic  collapse  and  densification  of  the 
foam  starts  at  the  impact  side  and  is  highly  localized,  while  the  rest  of  the  porous 
material  remains  elastic.  The  deformation  and  densification  originates  from  the 
impact  point  along  the  axis  of  penetration  and  propagates  in  the  direction 
perpendicular  and  traverse  to  the  applied  load.  Hence,  the  deformation  of  the 
aluminum  backing  plate  can  occur  only  when  the  foam  is  fully  compacted.  This 
delays  the  deformation  and  damage  growth  in  the  backing  plate  and  has  the 
potential  to  improve  ballistic  performance  of  the  armor. 

From  the  simulation  results,  the  predicted  behavior  of  the  laminated  target 
plate  under  ballistic  attack  is  similar  to  experimental  observations.  The  author 
observed  typical  brittle  fracture  failure  in  the  ceramic,  deformation  in  the 
Dyneema®  layer,  and  porous  compaction  of  the  PU  foam.  The  composite  plate 
resisted  complete  penetration  in  both  the  simulation  and  experimental  results. 
The  final  deformation  shape  and  deflection  measurements  from  the  simulation 
correlates  well  with  experimental  results  (see  Table  7).  Figure  44  shows  the 
experimental  results  for  Shot  13.  Ong  [2]  has  also  demonstrated  through  impact 
tests  that  the  ballistic  performance  of  the  Dyneema®  layer  degrades  severely 
without  the  high-strength  ceramic  layer.  Similarly,  the  ceramic  layer  alone  has 
poor  ballistic  performance  as  it  exhibits  extremely  brittle  fracturing  during  the 
impact  event.  Therefore,  it  is  necessary  to  adopt  a  multi-layered  armor  structure 
concept  for  improved  ballistic  performance. 


69 


PRESSURE  <kPa  )  PRESSURE  (kPa  ) 


Ceramic.  Al?0^ 


TIME  (ms) 


FR-6720  PU  Foam 


Figure  43.  Simulation  Results  of  Shot  Number  13 
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Figure  44.  Experimental  Results  of  Shot  Number  13 
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3. 


Ceramic-Dyneema®-FR6710-AI  (C-D-P1-AI)  Target  Plate 


For  this  system,  the  target  plate  consists  of  a  6  mm  ceramic  plate,  5  mm 
Dyneema®  layer,  and  5  mm  FR-6720  polyurethane  (PU)  foam  0.16  g/cm^  backed 
by  a  thin  Aluminium  6061 -T6  plate.  This  is  very  similar  to  the  test  above,  but  with 
a  lower  density  PU  foam.  The  measured  initial  impact  velocity  of  the  projectile 
was  481  m/s.  The  gauge  locations  in  the  simulation  setup  were  as  per  the  C-D- 
P2-AI  target  plate. 

In  general,  the  failure  mechanism  is  similar  to  the  C-D-P2-AI  target  plate. 
Compared  to  the  plate  configuration  with  the  more  dense  foam,  the  less  dense 
foam  more  strongly  attenuates  the  compression  shock  wave  (see  Figure  45). 
There  was  a  reduction  in  pressure  from  0.2GPa  at  the  last  gauge  point,  namely 
gauge  1 1  in  the  Dyneema®  layer,  to  65  MPa  at  gauge  16  in  the  PU  foam.  This  is 
equivalent  to  a  pressure  reduction  of  67.5%  across  5  mm  of  porous  foam. 
However,  the  author  observed  the  pressure  wave  rise  time  to  be  shorter  than  the 
denser  foam.  Hence,  the  damaging  compression  wave  reaches  the  backing 
plate  earlier  than  the  more  dense  foam  configuration,  and  backing  plate  damage 
deformation  and  damage  growth  occurs  earlier.  The  author  also  observed 
complete  crushing  of  the  foam  near  the  impact  point  and  he  observed  heavy 
erosion  in  the  simulations,  indicating  that  the  material  could  have  spalled  or 
pulverized.  This  probably  explains  the  higher  dynamic  deflection  obtained  for  the 
C-D-P1-AI  target  plate  deformation.  In  the  actual  porous  material,  micro-cracks 
can  initiate  in  the  brittle  porous  material  if  sufficiently  large  tensile  waves  are 
generated. 

Figure  46  shows  the  post-impact  condition  of  the  target  plate.  The  author 
observed  similar  failure  mechanisms  as  in  the  test  above.  This  composite  plate 
also  resisted  complete  penetration  but  with  a  larger  dynamic  deflection.  The  final 
deformation  shape  and  deflection  measurements  from  the  simulation  agree 
reasonably  well  with  experimental  results  (see  Table  7). 
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Figure  45.  Simulation  Results  of  Shot  Number  14 
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Figure  46.  Experimental  Results  of  Shot  Number  14 
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4. 


Ceramic-Dyneema®-Alulight®  AFS  (C-D-M)  Target  Plate 


The  author  modeled  the  target  plate  consisting  of  a  6  mm  ceramic,  5  mm 
Dyneema®’  and  12  mm-thick  Alulight  aluminum  foam  sandwich  (AFS).  He 
measured  the  projectile  as  having  an  initial  impact  velocity  of  465  m/s.  The 
gauge  locations  in  the  simulation  setup  were  as  per  the  C-D-P2-AI  target  plate. 

Similar  to  the  other  plate  configurations  with  porous  materials,  the 
simulation  results  showed  that  the  porous  aluminum  strongly  attenuates  the 
shock  pressure  (see  Figure  47).  There  was  a  reduction  in  pressure  from  0.2 
GPa  at  last  gauge  point,  namely  gauge  11  in  the  Dyneema®  layer,  to  28  MPa  at 
gauge  20  in  the  porous  aluminum  foam.  This  is  equivalent  to  a  pressure 
reduction  of  86%  across  5  mm  of  porous  aluminum  foam.  There  is  coupling  of 
the  strong  pressure  attenuation  with  widening  of  the  compression  shock  rise 
time,  which  then  introduces  significant  delay  in  wave  propagation.  The  air/gas 
filled  cellular  structure  of  the  aluminum  foam  makes  stress  wave  propagation 
difficult.  The  cell  wall  acts  as  a  wave  guide  and  disperses  the  stress  wave. 
Porous  compression  is  initiated  by  the  cell  wall  buckling  followed  by  plastic 
deformation  of  the  cell  wall,  which  then  leads  to  localized  densification.  The 
simulation  and  experimental  results  showed  that  the  plastic  collapse  and 
densification  of  the  foam  is  highly  localized,  while  the  rest  of  the  material  remains 
intact. 

Figure  48  shows  the  post-impact  condition  of  the  target  plate.  As  before, 
the  author  observes  similar  failure  mechanisms.  The  composite  plate  resisted 
complete  penetration  with  a  relatively  small  dynamic  deformation.  The  final 
deformation  shape  and  deflection  measurements  from  the  simulation  correlate 
well  with  experimental  results  (see  Table  7). 
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Figure  47.  Simulation  Results  of  Shot  Number  12 
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Figure  48.  Experimental  Results  of  Shot  Number  12 
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5. 


Comparison  between  Simulation  and  Experiment 


Table  7  is  a  summary  of  the  simulation  and  experimental. 


Overaii  Depth,  b 
(mm) 

Dynamic  Defiection 
(mm) 

Armor  Plate 
Configuration 

Plate 

Thickness 

(mm) 

Experiment 

Simulation 

Experiment 

Simulation 

C-D 

11 

45.97 

50 

34.97 

39 

C-D-P1-AI 

17.5 

35.12 

40 

17.62 

22.5 

C-D-P2-AI 

17.5 

32.64 

34.5 

15.14 

17 

C-D-M 

23 

33.57 

33.8 

10.57 

10.8 

Legend: 


Symboi 

Materiai 

Thickness  (mm) 

Ai 

AI  6061-T6  Aluminum 

1.5 

C 

AI  203-99.5  Ceramic 

6 

D 

Dyneema 

5 

M 

Alulight  aluminum  foam  sandwich  (AFS) 

12 

PI 

FR-6710  rigid  polyurethane  foam 

5 

P2 

FR-6720  rigid  polyurethane  foam 

5 

S 

AISI4140  steel 

4.76 

Table  7.  Comparison  of  Experimental  and  AUTODYN®  results 


The  AUTODYN®  simulations  produce  results  that  show  a  slightly  larger 
deformation  as  compared  with  experimental  results  due  to  the  confinement 
effects  of  the  experimental  setup.  In  the  actual  experimental  setup,  the  target 
plates  are  clamped  in  such  a  way  that  only  an  effective  area  of  50  mm  x  50  mm 
was  used.  For  the  author’s  simulations,  the  target  plate  is  simulated  as  100  mm 
X  100  mm  and  hence  produces  a  larger  effective  area.  Therefore,  the 
expectation  is  for  the  target  to  resist  penetration  for  higher  projectile  speeds  if  the 
total  effective  area  of  100  cm^  is  used.  However,  in  general,  the  characteristics 
of  the  simulation  results  agree  relatively  well  with  the  experimental  results. 
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From  the  above  results,  the  C-D  armor  plate  configuration  appears  to 
perform  worst  in  terms  of  dynamic  deflection  while  the  C-D-M  plate  configuration 
performs  best.  However,  the  ballistic  performance  cannot  be  ascertained  based 
on  the  experimental  results  alone  as  the  impact  tests  were  conducted  at  varying 
projectile  impact  velocities.  Since  the  above  simulation  results  agree  well  with 
the  experiment,  the  author  can  confidently  use  the  simulation  model  to  predict 
the  ballistic  performance  using  a  constant  projectile  impact  velocity  for  a 
comparison  study. 

6.  Ballistic  Performance  Comparison  Study 

In  order  to  perform  ballistic  performance  comparison  between  the  different 
armor  plate  configurations,  the  author  performed  AUTODYN®  simulations  using  a 
constant  projectile  impact  velocity  of  483  m/s.  Figure  49  shows  the  final 
deformation  shape  and  deflection  measurements  for  the  various  armor  plate 
configurations.  Figure  50  shows  the  time  taken  for  the  various  armor  plate 
configurations  to  detain  the  1”  cylindrical  steel  projectile  travelling  at  483  m/s. 
The  results  show  that  all  four  armor  designs  performed  credibly  with  no  projectile 
perforation.  All  four  armor  designs  performed  better  than  the  benchmarked 
AISI4140  steel  plate  (4.76  mm  thickness). 

From  a  ballistic  performance  point  of  view,  it  is  desirable  for  the  armor 
plate  structure  to  produce  the  least  dynamic  deflection  and  also  be  able  to  arrest 
the  projectile  at  the  shortest  possible  time.  From  the  simulation  results  (see 
Table  8),  the  C-D-M  armor  plate  configuration  has  the  best  overall  performance. 
However,  when  comparing  thicknesses  of  armor  designs,  it  is  important  to 
consider  the  dynamic  thickness,  in  other  words,  the  summation  of  static  and 
dynamic  thickness  [18],  which  is  equivalent  to  the  overall  depth  b.  The  areal 
density  or  mass  effectiveness  of  the  armor  design  is  also  a  critical  metric,  as  it 
will  affect  the  overall  weight  of  the  armor  design.  With  the  two  above 
considerations,  the  C-D-P2-AI  armor  plate  configuration  provides  the  best  overall 
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performance.  Incorporating  a  porous  layer  to  the  armor  design  improves  the 
performance  by  28%  compared  to  the  Ceramic-Dyneema®  (C-D)  composite 
armor. 


(a) 


(b) 


36  mm 


(c)  (d) 

Figure  49.  Simulation  Result  of  a  1”  Steel  Rod  Projectile  Impacting  (a)  C-D;  (b)  C- 
D-P1-AI;  (c)  C-D-P2-AI  and  (d)  C-D-M  Target  Plate  at  483  m/s 
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Figure  50.  Velocity  versus  Time  Plot  for  Projectile  Impacting  (a)  C-D;  (b)  C-D-P1- 

Al;  (c)  C-D-P2-AI  and  (d)  C-D-M  Target  Plate  at  483  m/s 
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Armor  Plate 
Configuration 

Areal 

Density 

(g/cm^) 

Overall 
Depth,  b 
(mm) 

Dynamic 

Deflection 

(mm) 

Time  taken  to 
arrest  projectile 
(ms) 

S 

3.737 

Complete  Penetration 

C-D 

2.771 

50 

39 

0.32 

C-D-P1-AI 

3.256 

40 

22.5 

0.20 

C-D-P2-AI 

3.336 

36 

18.5 

0.19 

C-D-M 

3.811 

34.5 

13 

0.20 

Equivalent  Areal  Density  to  C-D-P2-AI  Armor  Plate 

C-D-P1-AI 

3.336 

44.5 

22 

0.22 

Legend: 


Symbol 

Material 

Thickness 

(mm) 

Al 

Al  6061-T6  Aluminum 

1.5 

C 

Al  203-99.5  Ceramic 

6 

D 

Dyneema 

5 

M 

Alulight  aluminum  foam  sandwich  (AFS) 

12 

PI 

FR-6710  rigid  polyurethane  foam 

5 

P2 

FR-6720  rigid  polyurethane  foam 

5 

S 

AISI4140  steel 

4.76 

Table  8.  Ballistic  Performance  Comparison 


The  author  modeled  a  target  plate  consisting  of  a  6  mm  ceramic,  5  mm 
Dyneema®,  and  10mm  FR-6720  polyurethane  (PU)  foam  0.16  g/cm^  backed  by  a 
thin  Aluminum  6061-T6  plate.  This  armor  plate  configuration  has  the  same  areal 
density  as  the  C-D-P2-AI  armor  plate,  and  so  the  author  can  better  compare 
performance  of  the  two  foam  materials.  The  author  assigned  an  initial  impact 
velocity  of  483  m/s  to  the  projectile.  The  objective  of  this  simulation  run  is  to 
compare  the  ballistic  performance  of  two  porous  foams  with  equal  total  areal 
density.  Figure  51  shows  the  simulation  results,  and  it  clearly  shows  that  the  C- 
D-P2-AI  armor  plate  configuration  still  has  superior  ballistic  performance.  This 
result  is  unexpected  as  a  less  dense  foam  of  equivalent  areal  density  is  expected 
to  produce  larger  PV-work  due  to  larger  porous  volume  compression.  The 
simulation  results  showed  complete  crushing  and  severe  erosion  of  the  foam 
near  the  impact  point,  indicating  that  the  material  could  have  spalled  or 
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pulverized.  This  could  have  impeded  the  energy  absorption  performance  of  the 
porous  material  later.  There  needs  to  be  further  investigation  to  verify  the 
author’s  predictions  and  conclusions.  Microscopic  examination  of  the  impact 
damage  to  the  porous  foam  may  provide  valuable  insights. 


Projectile  Velocity  vs  Time 
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Figure  51 .  Simulation  Results  for  C-D-P1-AI  Armor  Plate  with  3.336g/cm^  Areal 

Density 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  current  study  has  shown  through  both  numerical  modeling  and 
ballistic  trials  that  a  multi-layered  armor  concept  consisting  of  an  initial  high 
strength  layer,  followed  by  a  second  wave-spreading  layer,  and  a  shock 
absorbing  third  layer  of  porous  material  has  significantly  improved  ballistic 
performance  and  energy  absorption  compared  to  conventional  steel  armor  of 
equivalent  areal  density. 

In  this  study,  the  author  developed  a  material  model  in  AUTODYN®  using 
the  porous  P-a  equation  of  state  compaction  model  to  describe  the  porous 
material’s  behavior  when  subjected  to  a  transient  load  such  as  projectile  impact. 
The  P-a  compaction  model  describes  initial  elastic  compression  of  the  porous 
material  due  to  elastic  cell  wall  buckling,  followed  by  plastic  deformation  of  the 
cell  wall  when  the  applied  pressure  exceeds  the  material  yield  strength.  This 
results  in  a  permanent  volume  change  to  the  compacted  porous  material.  If  the 
applied  pressure  continues  to  increase,  the  material  is  dynamically  compacted 
until  the  applied  pressure  reaches  the  fully  compacted  state,  which  coincides  with 
the  fully  dense  solid  material  Hugoniot.  Incorporating  the  porous  P-a  model  with 
the  other  material  models  developed  for  the  above  multi-layered  armor  concept, 
the  simulation  studies  conducted  have  shown  good  agreement  with  ballistic  trial 
results  in  terms  of  material  behavior  and  overall  ballistic  performance. 

The  role  of  the  high  strength  ceramic  first  layer  in  an  impact  event  is  to 

deform  the  projectile  and  strongly  dissipate  its  energy  before  it  reaches  the 

second  layer.  The  ceramic  layer  possesses  high  compressive  strength  but  is  not 

sufficient  to  stop  the  projectile  on  its  own  due  to  its  brittle  fracture  behavior.  The 

expectation  is  for  the  ceramic  material  to  cause  significant  plastic  deformation  of 

typical  small  arm  ball  ammunition,  which  is  made  of  relatively  soft  steel  core 

mass.  This  was  not  the  case  for  the  author’s  projectiles  since  they  were  made  of 
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relatively  hard  steel  that  did  not  deform  plastically  in  a  significant  manner.  The 
Dyneema®  wave-spreading  second  layer  utilizes  high  wave  speed  directional 
fibers  to  distribute  the  projectile  impact  energy  laterally  and  hence,  delays  the 
shock  wave  propagation  along  the  through-thickness  of  the  armor  material.  The 
third  layer  consisting  of  porous  foam  has  proven  to  be  a  good  shock  attenuator 
by  widening  the  shock  rise  time  to  delay  the  wave  propagation  and  strongly 
attenuating  the  shock  wave  by  absorbing  the  kinetic  energy  through  compaction 
of  the  porous  material  and  turning  it  into  waste  heat.  This  results  in  reduced 
dynamic  deflection  and  a  shorter  time  taken  to  arrest  the  projectile. 

The  current  study  has  also  clearly  proven  the  necessity  of  an  inertia 
backing  plate  to  confine  the  porous  material.  The  role  of  the  inertia  backing  plate 
is  to  prevent  the  porous  material  from  spalling,  which  impedes  its  ability  to  absorb 
the  impact  energy  efficiently.  Hence,  a  fourth  or  final  stopping  layer  consisting  of 
a  high  strength  or  penetration  resistant  material  will  be  essential  to  serve  the 
function  of  confining  the  porous  material  as  well  as  to  provide  the  final  attempt  at 
arresting  any  penetrator  that  has  penetrated  the  initial  three  layers. 

In  the  current  research,  the  author  investigated  three  types  of  porous 
materials.  Armor  solutions  containing  the  porous  material  layer  showed 
improved  ballistic  performance  compared  to  a  composite  plate  of  Ceramic- 
Dyneema®.  The  armor  plate  configuration  that  used  the  more  dense 
polyurethane  foam  had  the  best  overall  performance,  achieving  a  ballistic 
performance  improvement  of  28%  in  terms  of  dynamic  thickness  compared  to 
Ceramic-Dyneema®  composite  armor.  The  more  dense  foam  configuration  (C-D- 
P2-AI)  also  had  the  shortest  time  to  arrest  the  projectile.  Hence,  the  author  can 
confidently  conclude  that  the  attributes  of  porous  materials  contribute  positively 
to  the  weight/performance  ratio  of  the  armor  system  and  is  an  important 
constituent  of  the  multi-layered  armor  concept. 

The  current  study  has  also  revealed  through  analysis  of  the  simulation 
results  that  the  material  characteristics  of  each  layer  have  been  configured  to 

effectively  disrupt  and  absorb  the  kinetic  energy  of  an  incoming  penetrator. 
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Together  with  ballistic  trial  results,  the  author  concludes  that  the  sequence  of  the 
armor-layering  concept  is  fundamentally  correct  and  has  its  merits.  The 
enhanced  performance  with  weight  minimized  is  achieved  through  creative  role 
sharing  of  multi-functional  materials.  The  superior  performance  of  this 
innovative,  porous  foam  in  an  integral  layered  armor  system  is  a  step  forward  in 
the  design  of  more  damage-tolerant,  lightweight  personnel  and  vehicle  armor  for 
the  future  battlefield. 

B.  RECOMMENDATIONS  FOR  FUTURE  WORK 

Since  the  multi-layered  armor  concept  has  proven  its  benefits,  there  needs 
to  be  more  work  done  to  optimize  the  thickness  and  performance  of  each  layer. 
Obtaining  a  better  understanding  of  fundamental  stress  wave  propagation  in 
dissimilar  materials  and  material  failure  mechanisms  is  necessary. 
Recommended  future  work  includes  the  following: 

1.  Measured  Properties  of  Porous  Material 

Some  of  the  material  input  properties  for  simulation  input  were  derived 
empirically  using  limited  available  material  information  or  based  on  literature 
research.  It  is  important  to  perform  physical  experiments  to  derive  critical 
material  properties  such  as  bulk  sound  speed,  maximum  elastic  pressure  (Pe), 
and  solid  compaction  pressure  (Ps).  Ultrasonic  tests  can  measure  the  bulk 
sound  and  the  deduction  of  material  strength  properties  can  come  from  stress- 
strain  behavior  tests  conducted  at  high  strain  rates.  With  accurate  information  on 
material  properties,  it  will  be  possible  to  improve  the  fidelity  of  the  results,  and 
simulations  done  for  weight/performance  optimization  for  armor  designs. 

2.  Performance  Against  Ogive-shaped  Projectiie 

The  basis  for  current  simulation  models  and  experiments  are  a  flat-faced 
cylindrical  projectile.  One  should  perform  Numerical  modeling  using  a  typical 
ogive-shaped  projectile  to  replicate  actual  field  applications.  The  projectile  shape 
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is  likely  to  affect  the  material  behavior  in  an  impact  event  and  hence  the  overall 
ballistic  performance.  A  live  firing  experiment  could  validate  the  simulation.  The 
expectation  is  that  the  impact  velocity  of  a  typical  7.62  mm  ball  cartridge  will 
reach  as  high  as  900  m/s.  The  current  study,  due  to  limitations  of  the 
experimental  setup,  is  limited  to  impact  speeds  of  480  m/s  for  a  projectile  mass 
of  8.68  gm. 

3.  Effect  of  Pore  Structure  and  Pore  Size 

In  the  current  research,  porous  material  selection  is  limited  to  what  is 
commercially  available.  Based  on  the  material  characterization  results 
performed  on  the  polyurethane  foams,  the  two  different  foam  densities  possess 
similar  average  pore  size  and  pore  size  distribution.  Future  research  work  on 
porous  material  could  include  studying  the  materials  of  same  density  but  with 
different  pore  size  or  cell  structure.  The  study  will  be  able  to  provide  valuable 
insights  on  the  effect  of  pore  size  and  structure  on  energy  absorption  properties. 
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